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Abstract 

In this paper we prove a two-term asymptotic formula for for the spectral counting function 
for a 2D magnetic Schrodinger operator on a domain (with Dirichlet boundary conditions) 
in a semiclassical limit and with strong magnetic field. By scaling, this is equivalent to a 
thermodynamic limit of a 2D Fermi gas submitted to a constant external magnetic field. 

The original motivation comes from a paper by H. Kunz in which he studied, among 
other things, the boundary correction for the grand-canonical pressure and density of such a 
Fermi gas. Our main theorem yields a rigorous proof of the formulas announced by Kunz. 
Moreover, the same theorem provides several other results on the integrated density of states 
for operators of the type (— i/iV — /J-A.)^ in L^{Q) with Dirichlet boundary conditions. 
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C Existence and continuity of certain integral kernels |37 

1 Introduction and the main results 

1.1 The setting and the semi-classical results 

For h > we consider the self-adjoint operator 

Lh^{-ihV-Ay inL^ifl), 

where r2 C is a smooth, bounded and connected domain, and A is a smooth, real vector field 
defined on fi. Its curl, B := curl A, describes a magnetic field (which in the case of a simply- 
connected determines A uniquely up to the addition of a gradient field) . Throughout we shall 
assume that 

infi?(x)>0. (1.1) 

Our goal is to study the behavior of the trace Ti f{h~^Lh) in the limit ft — > for any Schwartz 
function /. Our main result will be a two-term asymptotic formula 



Tr /(ft-iL,0 = ft-i (Co(/) + h'/^C^{f) + o{h'^')) 



asft-^0 (1.2) 



with explicit coefficients Co(/) and Ci(/). The terms Co(/) and Ci(/) are bulk and boundary 
terms, respectively, and will be defined below. 

Due to the special structure of Co(/), asymptotics (|1.2p will imply two-term asymptotics for 
the number of eigenvalues of less than Eh, provided the energy E satisfies a certain gap 
assumption. In the particular case of a constant magnetic field, this assumption is satisfied for 
all energies different from the Landau levels. Thus we are looking for a magnetic analog of the 
formula [29l 



N{Eh,~h^A) ^ {Airh)-^ (\n\E -h^^'^\dn\E^/^ +o{h^^^)) as ft ^ , (1.3) 



where —A is the Dirichlet Laplacian on il and, for any self- adjoint operator A with discrete 
spectrum, N{X, A) denotes the number of eigenvalues less than A, counting multiplicities. We 
emphasize that, in contrast to our result in the magnetic case, (jl.3p requires some global geometric 
conditions on fl. Moreover, it is remarkable that our proof of a two-term asymptotics of the 
counting function is technically much less involved than the proof of (|1.3p . 

We emphasize that in (|1.2p we consider /(ft^^L/j) instead of f{Lfi), that is, we look at eigen- 
values of order ft instead of order 1. (For semi-classical results about Ti f{Ph) and N{E, Ph) for 
a large class of operators Ph we refer to [ISJ [TS] . These operators Pu are defined on the whole 
space, however, so that there are no boundary efi'ects.) The study of Tr f{h~^Lh) has applications 
to mean-field and thermodynamic limits in physical problems. It appears in |31| in connection 
with the mathematical justification of mean field theories for a large number of electrons in strong 
magnetic fields. Another motivation for us comes from the paper [30] which studies a Fermi 
gas confined to a large convex set submitted to a constant magnetic field. Boundary corrections 
are derived for the grand-canonical pressure, magnetization and density in the thermodynamic 
limit. Our main result (|1.2p provides a rigorous proof of this derivation and removes the convexity 
assumption. 
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We next describe the coefficients entering in (|1.2p . The term Co(/) is a bulk term, in the sense 
that it is of the form 

Co(/) := (2^)-i / fcfe(S(x), /)B(x) dx (1.4) 

involving the integral over n of the density 

bk{BJ):=f{{2k-l)B). (1.5) 

We emphasize that the leading order asymptotics Tr/(/i~^L/j) ^ h~^Co{f) are essentially known. 
Indeed, related formulas have appeared in [SJ [35] in the case of magnetic fields growing at infinity 
and in (311 [351 [Ml [3Z] in the presence of decaying electric potentials. 

Our main result concerns the second term. This is a surface term of the form 

Ci(/):=(2^)-i^ / Sfc(B(x),/)/BWrf'T(x), (1.6) 

depending only on the restriction of B to the boundary. Here da denotes integration with respect 
to the surface measure. In order to define the density Sk{B,f), we recall that the harmonic 
oscillator 

has eigenvalues 2fc— 1, fc = l,2,..., with corresponding normalized eigenfunctions 

Mt):-^;^{lM'^H,it)e^^, (1.7) 

where Hk is the fc-th Hermitc polynomial. For ^ G K we consider the Dirichlet realization of the 
operator 

^^ + (^ + *)' ini'(K+). (1.8) 

This operator has discrete spectrum with simple eigenvalues efe(^), fc = 1, 2, Let V'fc(-, denote 

corresponding normalized eigenfunctions chosen to be positive in a neighborhood of +oo. (This is 
possible since by the Sturm oscillation theorem eigenfunctions have only a finite number of zeros.) 
We put 

Sk{B, /) := r ( [ fiBckimMt, 0\'d^ - fiB{2k - 1))) dt . (1.9) 



Here and in the following, by a Schwartz function on [0, oo) we mean the restriction of a Schwartz 
function on R to [0,00). We are now in position to state our main result. The coefficient Ci(/) is 
well-defined according to the following 

Lemma 1.1. Let f be a Schwartz function on [0, 00). For B > the sum X^fcLi Sk{B, /) converges 
absolutely and is continuously dijferentiable with respect to B > 0. Moreover, for any = 1, 2, . . . , 

Sk{B,f)= [ (f{Bekm-.f{B{2k-l)) H \Mt)\' dt) d^ . (1.10) 



And here is the main theorem: 

Theorem 1.2. Let C be a bounded and connected domain with a boundary 951 assumed to 
be a finite union of disjoint, regular, simple and closed -curves. Let A G C^(ri,M^) and assume 
that B :~ curl A satisfies (|1.1[) . Then the asymptotics (|1.2p are valid for any Schwartz function f 
on [0,00) and with Coif) '^^'^ C'i(/) defined by (jl.4p . (|1.5p and (jl.6p . (jl.9p . respectively. 

In order to state the promised result about the eigenvalue counting function, we recall that 
A^(A, Lji) denotes the number of eigenvalues of Lh less than A, counting multiplicities. 
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Corollary 1.3. Let Q and A be as in Theorem M.'A and assume that for some positive integer K 



where B^i^ infxef2 -B(x) and S^ax := sup^gj^ B{-k). Then, for any E e ((2A' - l)i?max, (2A' + 

1)-Bmin), 

N{Eh, Lh) (k [ B(x) dx (1.11) 

^ f poo 

j,^]^ Jen J a 



. J{e B(x)e,(5)<-E} / 



as h ^ 0. In particular, z/ curl A = B is constant, then (jl.lip holds for any E \ {B(2k — 1) : 
A: = l,2,...}. 

We emphasize again that the crucial point in Theorem 11.21 and Corollary 11.31 is the existence 
of a second term in the asymptotic expansion. It is remarkable that a second term exists without 
any assumption on periodic orbits, in contrast to the case i? = [29]. The fact that this second 
term is a boundary contribution involving a one-dimensional auxiliary problem is a typical feature 
of two-term asymptotics; see, e.g., [34]. In [30] it is shown, using the alternative form of Sk{B, /) 
from Lemma 1 1.1 1 and the behavior of zeroes of special functions, that for any B > 0, the sum 

K 



fc=i 



1-/ iMt.ordn dt 



is a non-negative, decreasing function of i? G {{2K — 1)B, {2K + 1)B). It diverges to +oo as 
E \ B{2K — 1) and vanishes as E B{2K + \). This means that asymptotically the effect of the 
Dirichlet boundary is to reduce the number of eigenvalues, as compared to a bulk prediction. This 
is, of course, consistent with the non-asymptotic result from [20], which states that N{Eh,Lh) < 
{2TTh)-'^KB\n\ for any E e {{2K - l)B, {2K + l)B] and any /i > 0, provided the magnetic field 
B is constant and the domain Q, is tiling. Whether this extends to more general B and VL is an 
open problem. The corresponding problem for = is the famous Polya conjecture. 

Results similar to Theorem 11.21 and Corollarv ll.3l are derived in [191 [18] in the case of Neumann 
instead of Dirichlet boundary conditions and for / supported in [0,5] (or E < B \w (jl.lip ). In 
that case, however, Co{f) = 0, which makes the analysis considerably easier. 

A related interesting problem is the one in which the comparison operator lives on a strip 
of finite width [3J |3]. Enlarging the discussion, there exist some works dealing with spectral 
asymptotics for magnetic operators defined in the exterior of bounded domains, both with Dirichlet 
and Neumann boundary conditions (see |331 132j and references therein) . For a heuristic approach 
to sharp asymptotics for magnetic systems we refer the reader to |28j . 



1.2 The thermodynamic hmit 

In the case of a homogeneous magnetic field, the semi-classical limit is equivalent to the thermo- 
dynamic limit. We will now show how to recover and extend the results of |30| as corollaries of 
Theorem 11.21 For a set ft as in the previous subsection we define 

rih := h-^^^ O = {x e : /i^/^x £ fl} . (1.12) 

We fix a constant B > 0, put Ao(x) := ^{—X2, xi) and consider for any L > the operator 

:= {-iV - Aq)^ mL^{nL-2) 

with Dirichlet boundary conditions. Note that via the scaling x i-> L^^x, the operator i/^ is 
unitary equivalent to h'^Lh in L^(57) with h = L^^. Hence Theorem 11.21 implies 
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Corollary 1.4. Let fl be as in Theorem ] LSI let Ao(x) = ^{—X2,xi) for some B > Q and let f 
be a Schwartz function on [0, cxj). Then 

-L^ Trf{H^)^^f^b,{BJ) + L-^ (1.13) 

as L — > +00. 

Note that this coroUary is, in particular, vahd for f{E) = (e^(^-^) + l) and f{E) = hi(l + 
g-^i(£;-A')) for /? > and /i € R. These choices correspond to the Fermi-Dirac distribution and the 
grand-canonical pressure at chemical potential ^ and inverse temperature /3, respectively, which 
were considered in [5D] . In view of (jl.lOp , Corollary 11.41 provides a proof the formulas derived in 

1.3 Outline of the paper 

Our method to deal with the boundary contribution is reminiscent of the techniques developed in 
[23t [24l [TBI ITT] in the mathematical study of surface superconductivity. 

As in many other problems dealing with extended magnetic fields, gauge invariance and mag- 
netic perturbation theory plays a crucial role (see [7]-[l3], [2l]-[22] and references therein). 

Here is a short description of the content of this paper: 

• The proof of our main result. Theorem 11.21 is given in Section [SJ The main idea is to write 
the trace involving the original operator as a sum of traces involving model operators plus 
an error of order h°° . These model operators are either defined on the whole plane, or on 
slightly perturbed half-planes. The /i-asymptotics for the model operators are given in the 
previous three sections as independent results. 

• In Section [2] we consider operators defined on the whole plane. The main result here is 
Proposition [231 in which we give a sharp /i-expansion of the diagonal of the integral kernel 
oif{H). 

• In Sections[n|and[3|we treat the local trace ft,-asymptotics for operators defined with Dirichlet 
boundary conditions on the half-plane and on a perturbed half-plane, respectively, and when 
we restrict them near the boundary. The main technical results are Corollaries 13.21 and 14.21 
which, unlike Proposition [STTl are not formulated for integral kernels but for local traces. It 
would be interesting to prove a variant of these corollaries in terms of integral kernels. 

• Section [6| provides the proof of our results on the integrated density of states when the extra 
gap-condition is satisfied. 

• In the end we have three appendices containing various more or less known results, adapted 
and formulated according to our needs. 

2 Local asymptotics for operators on the whole plane 

In this section we consider functions of the operator (— iV — A)^ on L^(IR^). From Appendix [Cl 
we know that /((— iV — A)^) has a continuous integral kernel if / is a Schwartz function. In 
this section we provide an approximation for the diagonal value of kernel of /((— iV — A)^) which 
depends only on the value of the magnetic field at the point where the kernel is evaluated. This 
approximation is good provided the magnetic field is slowly varying. The precise statement is the 
following. 
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Proposition 2.1. Given f e 5(R), M > and N > 0, there is a CMnstant C > such that, for 
any B G (^^(M^) satisfying \\B\\^ + \\VB\\oo < M, for any xq G , h G (0, 1] and I G [1, oo) such 
that 

h-^/^ sup |VB(x)|+/i"i sup llHcss B(x)|| < M (2.1) 

|x-xo|<<? |x-xo|<£ 

and for any A, Aq G C'^(R^;M^) satisfying curl A = B and curlAg = -B(xo) one has 

|/((-zV- A)2)(xo,xo)-/((-2V- Ao)2)(xo,xo)| < C(/i + r~). (2.2) 

Later, in our proof of Theorem II. 21 we will apply this proposition with I = h^'^ for some e > 0. 
Then, when A'' is chosen large enough, we obtain an error of the order 0{h). The crucial point for 
our proof of Theorem 11.21 is that this is o(/i^/^). We believe that the order 0{h) is best possible. 
Note also that we do not require the magnetic field B to be positive. 

Before beginning with the proof, we will derive from the pointwise estimate in (|2.2p a weaker, 
integral estimate for Tr g/((— iV — A)^) with g of 'small' support. To do so, we also recall that if 
curlAo = i? is a non-zero constant, then 

/((-zV - Ao)2)(xo, xo) = ^ E - = ^ E ^'^■d^l' /) (2.3) 

fc=i fe=i 

for all xo G M^. It is easy to see that the right side of the above formula can be continuously 
extended to _B = by setting it equal to (47r)~^ /(A) dX. 

Corollary 2.2. Given f G 5(M), M > and N > 0, there is a constant C > such that, for any 
B G C^(M^) satisfying ||i?|loo + 1| Vi?||oo < M, for any xq G R^, h G (0,1] and £ G [l,oo) satisfying 
([2J|) . for any g G L°°{{\ ■ -Xo| < £/2}) and for any A G C3(R2;M2) ^-^f^ curl A = B, one has 



oo „ 

Tr[g/((-zV-A)2)] / 5(x)6.(|S(x)|, /)|i3(x)| dx 



<C\\g\\^e' (h + (2.4) 



Indeed, this is an immediate consequence of (|2.3p and (|2.2p noting that (12.11) holds for any 
X G supp g with £ replaced by i/2. 

Wc now turn to the proof of Proposition 12.11 We shall use the functional calculus based on 
the Helffer-Sjostrand formula, which we recall next. Assume that / G ^(M) is real-valued. For 
any given iV > 1, we can construct an almost analytic extension fa.N which obeys the following 
conditions pS] : 

1. fa.Niz) = /(z), Vz G M; 

2. supp(/a,Ar) C {z G C : |Imz| < 1} =: V; 

3. fa.N e S^(V); 

4. there exists Cn > such that for all z G I? 



dfa.N , . 



llmzl^ 



dz 

Let a be any self-adjoint operator. Assuming that N >2, the Helffer-Sjostrand formula reads: 
f{H)^- f ^^{H~z)-Uxdy, with z = x + iy. (2.6) 

TT J J, dz 

The proof of the following lemma is elementary. 
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Lemma 2.3. Assume that B e C^{W) satisfies \\B\\ao + ||VB|loo < M and 

sup |VB(x)| + h^^ sup II Hess B(x)|| < M 

|x|<f |x|<£ 

for some positive M,h,£. We define 

A(x):=^ dttB{tx){~X2,Xi), Ao(x) :^(-X2,Xi), 



(2.7) 



and (5A := A — Aq. Then for all x G 



p2 



|(5A(x)| < A/|x| and | div (5A(x)| < iAf|x| . (2.8) 

3 



Moreover, for all x with |x| < ^ 



|(5A(x)| < iM/ii/2|x|2 |div5A(x)| < ^M/i^/^i^l (2.9) 



<5A(x)-i(VB(0)-x) (-2:2, a;i) 



<iA//i|xp. (2.10) 



Proof. The first bound in p.Sp follows by writing 

6A{x)^[ dtt{B{tx) - B{0)){~X2,Xi) 
Jo 

and estimating |i?(tx) — -B(0)| < 2A/. Similarly, for the second bound we write 



(2.11) 



div (5A(x) = / dtt^ W B{tyi) ■ {-X2,xi) (2.12) 
Jo 



and estimate |Vi3(tx)| < M. 

For the first bound in p.9p we use that B(tx) - 5(0) = J^ds\/B{ sx) • X in (f2lT|) to obtain 



1 

^A(X) = - y ds (1 - [VS(SX) • X] {~X2, Xi) 



(2.13) 



The claimed inequality then follows from |Vi3(sx)| < Mh^/"^. For the second bound we use the 
same estimate in p.l2p . 

Finally, by ([2T3l) we have 

5A(x) - i [VB(0) • x] (-X2, xi) = i ^ ds (1 - s2) [(VB(sx) - VB(0)) • x] (-a;2, xi) 

and, since VB{sx) — VB{0) = da Hess B{ax) x, 

1 1 /"^ 2 \ 

(5A(x) - - [VB(0) -x] (-X2,.Ti) = - y ^'^ (^^ " + 3^^^ j (x-HessB((7x)x)(-a;2,xi). 

In view of ||Hess_B((Tx)|| < Mh, we arrive at (|2.10p . □ 



Proof of Proposition \2.1\ Step 1. Representation via the Helffer-Sjostrand formula 

We may assume that xq = 0. Moreover, since the diagonal of the integral kernels are gauge- 
invariant, we may assume that A and Ag are given by (|2.7p . We write H := (— iV — A)^ and 
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Ho := (— iV — Aq)^ in L^(R^). Applying the resolvent identity twice, we find that for any 
^ G C \ [0, oo) (the equality holds on functions with compact support): 

{H - z)-^ = {Ho - z)-^ - {Ho - z)-^W{Ho - z)-' + {Ho - z)'^W{H - z)-^W{Ho - z)-^ 

where 

W -{-iV - Ao) -SA-SA- {-iV - Aq) + \5A\^ 

and 5A := A — Aq. 

The above identity needs some justification since W is unbounded, and this is what we 
do next. Let us consider one piece building the product W {Hq — z)~^ , i.e. the operator 
5 A ■ (— iV — Aq) {Ho — z)"^, and let us apply it on a compactly supported function Then we 
can write: 

5A ■ {-iV - Ao) {Ho - z)-^^Jj =e-^l-l{e-|7l-l5A} 

• {e^l-lHV - Ao) {Ho - z)-ie-^l-l}{e^l--"l^}. (2.14) 

Using the general estimates established in (|B.13|) and (jB.14[) it follows that the resolvent sand- 
wiched with exponentials is bounded, and we also gain some exponential decay to the left. Half of 
it is used to bound the linear growth in SA, while the other half can be used as a priori decay in 
case we have other W^s and resolvents to bound. Moreover, the norms grow at most polynomially 
in |Re2;| and in l/|Imz|. 

Hence by the Helffer-Sjostrand formula for any almost analytic extension fa of / satisfying 
the decay condition (|2.5p . we can write the following operator identity which holds at least on 
compactly supported functions: 

f{H) - f{Ho) = -- I ^{Ho- z)~^W{Ho - z)"i dxdy (2.15) 

TT J-p dZ 

+ - [ ^{Ho~z)-'W{H~z)-^W{Ho~z)-^dxdy 

TT Jjj dz 

where T) = {z G C : |Im zj < 1}. Hence we need to prove estimates on the diagonal of the integral 
kernels of {Ho - z)-'^W{Ho - z)-^ and {Ho - z)~^W{H - z)-^W{Ho - z)-'^ . In order to do this, 
we need use the estimates from Lemma 12.31 on 5 A and its divergence. In the remainder of this 
proof we always assume that z Cz V. 

Step 2. The term {Ho - z)-'^W{Ho - z)-^ . We write 

{Ho - z)-^W{Ho - z)-i = - {{-iV - Ao){Ho - zY^ ■ SA{Ho - z)-' 

~ {Ho - z)-HA ■ (-zV - Ao){Ho - z)-^ 
+ {Ho ~ zr'\5A\\Ho - zr\ 

Since the first term on the right side coincides with the adjoint of the second term after replacing 
z by z, it suffices to consider the second and the third term. 

We begin with the second term. Let Qo{-,-,z) be the integral kernel of {Ho ~ z)~^ and Qi{-,-,z) 
be the (C^-valued) integral kernel of (— iV — Ao){Ho — z)~^. Then the integral kernel of {Ho — 
z)-^SA ■ (-iV - Ao){Ho - z)-^ is 

/ dy5o(x,y,z)<5A(y).gi(y,x',z). (2.16) 

Since 

|go(x, x', z)\ < C{z){l + I In |x - x'lDe-l'^-'^'l/^^^' , |gi(x, x', z)\ < C(z)|x - x'p 'e-l^-^'l/^^^^ 
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by Propositions IB.8] and IbTqI and since |5A(x)| < Af|x| by Lemma [2.31 this integral converges 
absolutely for any z G C\M and for any x and x'. Moreover, one can show that the above integral 
depends continuously on (x, x'). 

Therefore we can and will set x = x' = in the following. In order to estimate the integral 
(|2.16p we split it according to whether |y| < £ or not. In the first case we note that Qo{0,y,z) 
depends on |y| only. This follows from Mehler's Formula (see (|B.2ip '). Hence the tangential 
component of CJi(y, 0,z) is also a radial function times (— 7/2,2/1)- This implies that for any fixed 
vector c G M^, 

go(0,y,z) (c-y) (-ya, 2/i) • ^iCy, 0, z) (2.17) 
is odd with respect to y — y, and therefore its integral over {|y| < £} is zero. Hence by (|2.10p . 



dyeo(0,y,z)(5A(y) •gi(y,0,z) 



|y|<^ 



< hlh [ dy 1^0(0, y,z)||ynei(y,0,z)| 
8 J\y\<e 



By Propositions IB. 81 and IB. 91 there are C and S such that for all z G I? the latter integral is 
bounded by 

C— / dy (l + |ln|y||)|y|2e-^lyl. 



Recall that we are using the notations r = a/ (Re z) and 77 = |Imz|. One easily sees that for any 
e > there is a Ce such that the latter is bounded for all z GV by 



^10 / J, ^ 



In the region where |y| > £, we use the first estimate in (|2.8p to get 

/ dy GoiO^y, z)S Aiy) ■ g,{y,0, z) <M f dx |^o(0, y, z)| |y| |^i(y, 0, z)| 
J\y\>e J\y\>e 



(2.18) 



By Propositions IB.8] and |BT9| there are C and S such that for all z with |Imz| < 1 the latter 
integral is bounded by 



^10 r 

C— rfy (l + llnlylDe-^l^'l. 

r J\y\>i 

One easily deduces that there is a C such that the latter is bounded for all z G I? by 

(This can certainly be improved.) Using t^e~* < (N/e)^ , we conclude that for any > there 
is a Cn such that for all z G the latter is bounded by 

Cn—(-) 

This, together with (|2.18p . proves that for any iV > there are Cn and L such that for all z E T), 

L 



f dygo(0,y,z)<5A(y).gi(y,0,z) 

JR2 



<Cn{-] + 



Next, we briefly turn to the term {Hq — z) ^\6A\'^{Hq — z) As before one shows that this 
has a continuous integral kernel and that for any iV > there are C^v and L such that for all 
z G X» 



dyeo(0,y,z)|5A(y)|25o(y,0,z) 



< 
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This proof uses the first inequalities in (j2.8p and (|2.9p . The latter leads immediately to a power 
(/i^/^)^ = h and therefore one docs not need a symmetry argument in this case. 

Since for any given L, the almost analytic extension fa in (j2.15p can be chosen to satisfy for 
some constant Cl, 



\ L 



az 

for all z G V, the above estimates prove that the integral kernel of 

^(Ho - z)-^W{Ho - z)-i(0,0) dxdy 

at (0, 0) exists and is bounded in absolute value by Cn {h + £^^) for any N. 
Step 3. The term [Hq - z)-'^W {H - z)'^W {H^ - z)'^ . 

Using that (— iV — Aq) • = 5 A ■ (— iV — Aq) — i divJA we have two representations 

W = -2{-iV - Ao) -SA + V* 

and 

W = -25 A ■ (-iV -Aa) + V 
where V := |'5Ap + idiv 5 A. This allows us to write 

{Hn- z)-^W{H - z)~^W{Hn- z)"^ (2.19) 
=4 ((-^V - Ao)(ffo - z)-^)* ■ SA{H - z)-^6A ■ (-zV - Ao)(i?o - ^)"' 

- 2(i/o - z)-'V*{H - z)-'dA ■ i-iV - Ao){Ho - z)-' 

- 2 ((-*V - Ao)(i7o - zy'T ■ SAiH - z)-^V{H^ - z)'^ 

+ {Ho - z)-^V{H - z)-'V{Ho - 

We begin by considering the first summand. If Q{-, •, z) denotes the integral kernel of (H — z)^^ , 
then the integral kernel of ((-iV - Ao){Ha - • SA{H - z)-^6A - (-iV - Ao)(i?o - z)~^ is 

given by 

f dy f dy' g,{y, x, z) • SA{y) g{y, y' , z) 5A{y') ■ (y', x', z) 

As before, it follows from Lemmas IB. 81 and IB. 91 and the first estimate in (|2.8p that this integral 
converges absolutely and depends continuously on (x, x'). 

Moreover, the same argument as in Step 2 shows that the above integral, evaluated at x = 
x' = coincides with 



/ dy f dy'g,{y,0,z)-SA{y)g{y,y',z)SA{y')-g,{y\0,z) 
"'|y|<f •'\y'\<e 



up to an error which for any N can be bounded by Cat £ ^ for some Cn and some L = L{N). 

According to the first estimate in (|2.9p this integral is bounded by 

hi'h f dy f dy'|C?i(y,0,z)||yn(?(y,y',z)||yf |gi(y',0,z)|. 

According to Propositions IB. 81 and IB. 91 this double integral is bounded by 

(7^/ dy f dy'e-^lyl|y|(l + |ln|y-y'||)e"^|y-'''l|y'|e-^|y'l. 

/ X 4+£ 

One easily checks that for any e > this integral is bounded by I ) 



The other summands on the right side of (|2.19p are treated similarly. For the integral kernel 
of {Hq — z)^^V{H — z)^^V{Hq — z)^^ one again first localizes to |y'| < £ and |y| < i and then 
estimates V[y) by (|2.9p . The two factors of V lead to a factor of (/i^/^)^ = h. In the second and 
the third summand in (|2.19|) one can apply the Schwarz inequality and use the estimates for the 
first and the fourth summand. □ 
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3 Local asymptotics for operators on the half-plane 



In Corollarv l2.2l we studied Trg/((— zV — A)^) for the operator (— iV — A)^ defined on the whole 
space. Our next goal is to study similar traces for operators defined on a domain with boundary 
and for functions g supported close to the boundary. Once again, the idea is that if B varies slowly 
(measured in terms of a parameter h) on the support of the function g (measured in terms of a 
parameter £), then an operator with constant magnetic field provides a good approximation. In 
contrast to Corollarv l2.2[ however, the presence of a boundary leads to the appearance of a second 
term involving an integral over the boundary. As a first step towards proving this, we consider in 
this section the case of the half-plane = {x e : X2 > 0}- 

Proposition 3.1. Given f E S{M.), M > and N > 0, there is a constant C > such that for 
any B G C^{R\) satisfying ||B||oo + !| VB||oo < M, for any G Rx {0}, h e (0,1] and £ e [l,oo) 
such that 

sup \VB(x.)\ < Mh^/^ , (3.1) 

xSK^ , \x-xo\<l 

for any A G C'^{M.\,M?) and Aq G C^(R^,M^) satisfying curl A = B and curlAo = i?(xo) and 
for any g G i°°({| • — xo| < ^/2}) one has 

|Trg/((-*V - A)2) - Tr.g/((-zV - Ao)^)| < Cf\\g\\^{h^/^e^ + + r^) . (3.2) 

Here (-iV - A)^ and {-iV - Aq)^ are the Dirichlet realizations of the corresponding operators 

We emphasize that, in contrast to Proposition l2.H we only impose conditions on B and its first 
derivative. Therefore the result is not as precise as those from Proposition 12.11 and Corollarv 12.21 
This, however, will be sufficient for the application we have in mind. In our proof of Theorem 1 1.21 
we will choose £ = for some e > 0. If e < 1/6 and N is chosen large enough, then the error 
will be o{£) = o{h~'^), whereas we shall compute Tr g/((— iV — Aq)^) up to order 0{£). 

Before turning to the proof of Proposition 13 . 1 1 we deduce the following consequence. 

Corollary 3.2. Given f G 5(R) and M > 0, there is a constant C > with the following property: 
Let B G C^{M.l) satisfy \\B\\^ + ||VB||oo < M and let xo G R x {0}, h G (0,1] and £ G [l,oo) 
satisfy 

sup |VB(x)| < A//ii/2 and inf B{x) > M'^ (3.3) 

xGR^ , |x-xo|<C xSR^ , |x-xo|<^ 

Finally, assume that g G C^({| • -xo| < £/2}) satisfies 

sup \g{x)\+i-'^ sup |V5(x)|<M. (3.4) 

xSR^, |x-xo|<f/2 xGK^, |x-xo|<£/2 

Then for any A G C^(R^,R^) with curl A = B one has 

Trgf{{-iV-Af)-^Y. I 9i^)bk{B{^), f)B{^) dK 



-^ZJ / 0)sfc(B(xi, 0), /)v/S(a-i, 0) dxi 

^ k=i •'^ 

< C£^ {h}''^£^ + h£^ + r^) . (3.5) 

Here (— iV — A)^ is the Dirichlet realization of the corresponding operator in L^(R^). 

Remark 3.3. There is nothing special about the 1/2 in the support condition on g in Proposi- 
tion [3lT] and Corollarv 13.21 Indeed, the proof below works equally well when the support of g is 
contained in {| • — xo| < 9£} for a fixed (^-independent) G (0, 1). Of course, then the constant 
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will also depend on 6. Similarly, there is nothing special about the lower bound 1 on Again the 
proof works equally well when I S [^OiOo) for a fixed (/i-independent) > 0. The constant will 
also depend on t^. 

Proof of Corollary VJ.SX Our goal is to replace the term Tt gf{{~iW — Aq)^) in Proposition 13 . II by 
the two integrals on the left side of (|3.5p and to control the error. To do this, we first note that 
the operator (— iV — Aq)^ can be diagonalized explicitly. By gauge invariance we may assume that 
Ao(x) = b {—X2, 0) with b = B{x.q). Applying a Fourier transform with respect to the xi variable, 
we arrive at the operator —-^ + [pi + bx2)^- Via the change of variables t = b^X2 this operator 

is unitarily equivalent to b (^—^ + (&^'pi + t)^^ . This is the operator (|1.8p with ^ = b^^pi. We 

conclude that the eigenvalues and eigenfunctions of the operator — ^ + {pi + bx2)^ are given by 

bej(b^ipi) and biipj{b^X2,b~^pi). Hence 

Tr5/((-^V - Ao)^) = — E / , ffW^^C^ /) (3-6) 

+ ^ E / d^i dx2 ff(x) 1^ 6i/V(&efc(fe-^pi))|V^fe(&^X2, - bf{b{2k - . 

Since the function A ^ /(A(2fc — 1)) is in C^(0, oo) and since VB is controlled on the support 
of g by Mh^/"^. we have the estimate 



5(xK-(6,/)dx-^f; / 5(x)6,(i?(x),/)i?(x)dx 



< C h^'^l^ . 



(3.7) 



Here we also took into account that | suppgj < (7r/8)£^. 

The second term on the right side of (|3.6p needs more care. After rescaling X2 and pi it equals 



— V / dxi / dig(.Ti,tfe-5)fei/ 



/ fibekmMt,C)\'dC-fmk-l)) 



First, we replace g{xi, tb 2 ) by g(xi, 0). Because of the bound on the derivative of g we have 

\g{xi,tb-^)-g{xi,0)\ < Mi-Hb-i < M^/^rH. 
In Lemma I A. 3 1 we prove that 



E 

fc=i 



dtt 



fibekmMt,0\'dC~ f{{2k-l)b) 



< 00 . 



This allows us to replace (7 2) hy g{xi,0) and then to carry out the t integration. We obtain 

p POO f r 

J2 dxj dtg{xi,tb-i)b'^U / fibekmMt,0\^d^-f{{2k-l)b) 
j^^^ Jr Jo IJr j 

-E / dxig{x,,0)b'^hk{bj) 



< C. 



Finally, we use the fact that the function A 1-^ J^k f) is in C^(0, 00), see Lemma [A. 11 This, 
together with our assumptions on B and shows that 



00 00 

E / dxig{xi,Q)b^'^su{bJ)-Yl / 0) 0)sfc(B(xi, 0), /) 



(Note that in the last step we used i?(0) > 0.) This leads to the bound stated in the corollary. □ 



12 



We now turn to the proof of Proposition 13. II The following lemma is the analogue of Lemma 
13.41 We omit the proof, which follows the same lines. 

Lemma 3.4. Assume that B e C^(M^) satisfies ||i?||oo + ||Vi3||oo < and 

sup |VB(x)| < 

for some positive M, h, £. We define 

A(x) / dtB{xutx2){-X2,Q), Ao(x) 5(0) (-X2, 0) . (3.8) 
Jo 

and SA := A - Aq. Then for all x G 

|5A(x)| < 2Mx2 and \ div SA{x)\ < Mx2 . (3.9) 

Moreover, for all x G IR^ with |x| < £ 

|(5A(x)| < mi/2^2|x| and | div (5A(x)| < M/i^/^^.^ . (3.10) 

Proof of Proposition lS. l\ Step 1. Representation via the Heljfer-Sjdstrand formula 

We may assume that Xg = 0. Moreover, since the local traces are gauge-invariant, we may 
assume that A and Aq are given by p.Sp . We write H := (— iV — A)^ and Hq := (— iV — Aq)^ in 
L^(M^) with Dirichlet boundary conditions. With a justification as in the proof of Proposition l2.1l 
one has the resolvent identity in the form 

{H - z)-i = (Ho - z)-' - {Ho - zy'WiH - z)-' 

where 

W := -(-iV - Ao) -SA-SA- (-iV - Aq) + \5A\'^ 

and 5 A := A — Aq. Hence by the Helffer-Sjostrand formula for any almost analytic extension fa 
of / we can write: 

f{H)-f{Ho) = -- f ^iHo-z)-'WiH-z)-Uxdy, 

TT Jj, CfZ 

where again 23 = {z e C : |Imz| < 1}. In contrast to the proof of Proposition 12. 11 however, this 
identity is not sufficient for our purposes, since the operator {Hq — z)~^W{H — z)^^, even when 
multiplied by cut-off functions from the left and from the right, is not trace class. Our way out 
is to note that the above integral does not change if we subtract an analytic (operator-valued) 
function to the integrand. We subtract {Hq + 1)^^W{H + 1)^^ and compute 

{Hq - z)-^W{H ~ z)-^ - {Hq + \)-^W{H + l)-i 

= (z + 1) {{Hq - z)-\Hq + l)-^W{H - z)-i + {Hq + l)-^W{H + 1)-\H - z)-^) . 

Thus 

f{H) - f{HQ) =-- f dxdy ^{z + 1) 
Jv oz 

X {{Hq - z)-\Hq + l)-^W{H - z)-i + {Hq + l)-^W {H + l)-\H - z)-^) . 

Next, we denote by 77 the characteristic function of the disk {|x| < I}. Then, by cyclicity of 
the trace and the fact that rjg = g, we have 

TTg{f{H)~f{HQ)) = --f dxdy^{z + l) 

X Tr {i^{Hq - z)-\Hq + l)-^W{H - z)-^g + g{HQ + 1)" V(F + 1)-\H - z)-\) . 
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In the remainder of this proof we will show that the operators //(-ffo — z)^^{Ho + 1)~^W{H ~ z)~^g 
and g{HQ + l)~^W{H + 1)~^{H — z)~^r] are trace class and we derive bounds on their trace norms. 
We begin by bounding 

UHo-z)-\Ho+ir'W{H~z)-'g\\s^^L.^^.^)) < UHo-z)-\Ho+l)-'\\B,(^LHMl))\\W{H-z)-'g\\ 
and 

\\g{Ho+l)-'W{H+l)-\H-z)-W\B,^LHMl)) < \\{H+l)-\H-z)-'7^\\B^^i^.^^2^^^\\g{Ho+l)-'W\\ . 
The trace class factors on the right side are controlled by 

h{Ho + ir^B,(L-(Rl))<Cf and 1| (i/ + 1)-2^||b,(l2(r2 )) < _ (3^^^) 

These bounds follows from Proposition IB. 51 together with the fact that one needs 0{£'^) squares 
of unit diameter to cover the support of 77. In addition, we use the fact that {Hq + l)(i?o — = 
1 + {z + I) {Ho — z)^^, which implies 

\\{Ho + l){Hn- z)-^\\ < l + \z + l\\lmz\-K (3.12) 

A similar bound holds for H. 

To summarize, we have shown that 



\TTg{f{H)-f{Ho))\<C£^ [ dxdy 



dfa 



cfz 



i. + iiii + " 



Im z 



X {\\W{H - z)-'g\\ + \\g{Ho + l)-'W\\) , (3.13) 

We now decompose W into a part close to the origin and a part far from the origin. More 
precisely, we write 

W = riW+il~ g)W = Wi] + W{1 - r;) . 

We employ the two different forms for the two different operators we have to estimate. In the 
following two steps of the proof we will derive the bounds 



r2 

■q 

and 



2 

\\t^W{H ~ z)-^g\\ + |l.g(Fo + l)-^Wii\\ < C\\g\\Jy {h^'H^ + M^) (3.14) 



J.6+N 

11(1 - g)W{H - z)-'g\\ + \\g{Ho + 1)-'W{1 - r/)|| < CMo. ^ (3.15) 



for z E v. We recall the notations r = \J (Re z) and 77 = |Imz|. The constants C and Cjv in p.l4p 
and (|3.15l) depend only on M and the latter one also on > 1. These bounds, when inserted 
into (|3.13p . are integrable with respect to z (provided the decay order A in (|2.5p of the almost 
analytic extension is large enough) and lead to the behavior claimed in the proposition. 

Stej) 2. Bounds on rjW and Wrj. The idea behind our estimate in this case is that the 
coefficients in the operators gW and Wg are non-zero only close to the origin, and even there, 
they are small, as quantified by Lemma [3.41 We shall show that 



\\{Ha + l)-^/^Wg\\oo<c(^h'/^e^ + M'') and \\gW{H + < C (^h^/^i^ + M'') (3.16) 

with a constant C depending only on AI. This, together with p.l2p . implies p.l4p . 
To prove (j3.16p . we compute 

Wg = -2{-idi - Ao,i)gSAi+gV* , 
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where again V* ~ (i^^i)^ ~ idiSAi. Here Aq.i denotes the first component of Aq and similarly 
for SAi. The Schwarz inequality implies that 

Lemma [3.41 together with the properties of 77 implies that ||77(5Ai||oo < 4A//i^/^£^ and ||77V^||oo < 
16A/2/i£4 2Mh^/^e. Thus 

Wo^p)\ < c [h'^'e' + he^) {Ho + iW^'m' 

with C = 8M + 16AP + 2M. This proves the first inequality in p.l6p . The second one follows in 
the same way by using the representation 

rjW = -2rjSAi {^idi ~ Ai) ~ r/V* . (3.17) 



Step 3. Bounds on (1 — 77)!^ and W{1 ^ if). The idea behind our estimate in this case is that 
the coefficients in the operators (1 — ri)W and W{1 — rj) vanish on the support of g. This fact, 
combined with the exponential decay of the resolvents, leads to errors which decay faster than any 
polynomial in l. To make this precise we apply Proposition IB. 71 to the operators (1 — rj)W and 
W{1 — rj) and deduce that given M there are constants C and (5 > such that for any z £ V and 
for any characteristic functions x a-nd x of disjoint sets of diameter < 1 and of distance d > 1 one 
has 



||X(1 - V)W{H - z)-^x\\ < C-e-^'"'/nix(l - v) 

V 



The fact that C only depends on M follows from the fact that the left side of (|B.13p can be 
bounded in terms of M and ||x(l ^ v)^'^\\oo- Indeed, by the argument following (jB.13[) it suffices 
to bound the norm of e"^'~^°'^ xi^ ~ 'r])We^"^'^^"'^ {Hd — in terms of these quantities. This 
follows easily from the representation (cf. p.l7p ) 

(1 - 7j)W = -2(1 - Tj) SAi {~idi - Ai) - (1 - i-i)V* , (3.18) 



together with the bounds from Lemma YdM on 5Ai and V . 

We now choose the Xj 's to be characteristic functions of sets of diameter < 1 such that Xj 
is the characteristic function of the support of 1 — 1]. (For instance, let Xj be the characteristic 
function of a square on the lattice (2~^/^Z)^ times the characteristic function of the support of 
1 — 77.) Similarly, let the Xfe's be characteristic functions of sets of diameter < 1 such that Xk 
is the characteristic function of the support of g. Note that since g is supported in {|x| < £/2} 
and 1 — 77 is supported in {|x| > ^}, the distance dj^k between the supports of Xfc and Xj is at least 
Thus the above bound allows one to deduce (for £ > 2 at least; the remaining case is easily 
handled) 



\\{l-fj)W{H~z)-^g\\ 



1/2 / \ 1/2 



C\\g\\oo- suplix.x^lloo^e-*'"^-'-/'- sup^p-*"'^-^/'^! 

V ^ k / V '= J 

Using elementary estimates one can bound the product of the two square roots by a constant times 
(-{r /rffe~^^^/'^ for some Q <5 < 5. Given TV > we use the bound i^e~* < {N/e)^ to arrive at 
the first bound claimed in p.isp . The proof of the second one is similar. This concludes the proof 
of Proposition [3?n □ 
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4 Local asymptotics for operators on a perturbed half-plane 

We now extend the results from the previous section to 'perturbed half-planes'. By this we mean 
domains of the form 

V = {yeR^: V2>q{yi)}- (4.1) 

with a function q G C^(R). We note that the change of variables y i-> x = (2/1,2/2 — 9(2/1)) 
transforms V into M^. In this section, we will typically use a~ to denote an object on which 
corresponds to the un-tilded object on V . What 'corresponds' means depends on the object we 
are dealing with. For scalar functions g onV we define functions g on by 

5(x) = g{xi,X2 + q{xi)) . 

For vector fields A e C'^{'P,R'^) wc define A on by 

A(x) {Ai{xi,X2 + q{xi)) + q'{xi)A2{xi,X2 + q{xi)), A2{xi,X2 + q{xi))) . 

This definition is motivated by the fact that curl A(x) = curl A(a:i, X2+q{xi)), that is, A generates 
the same magnetic field as A in the new coordinates. Finally, ii H = (— iVy — A)^ is a magnetic 
operator in L^{V), with Dirichlet boundary conditions, say, then 

H=(-id^,-Ai- q' (-id.,,, - ^2) ) ' + i-idx, -A2^ 



denotes the corresponding operator on L^(K^). The operators H and H are unitarily equivalent 
via the above change of variables. 

The next proposition, which is the main result of this section, quantifies in which sense H is 
close to the simpler operator (— iV — A)^ when the boundary is almost flat. 

Proposition 4.1. Given f e 5(]R), M > and N > 0, there is a constant C > such that for 
any q G C^iR) satisfying \\q'\\^. + \\q"\\oo < M, for any B G C[R\) satisfying \\B\\^ < M, for 
any yo G dV , h G (0, 1] and I G [1, M/i^^/^J such that q'{ya,i) = and 

sup \q"{y^)\<Mh^/\ (4.2) 
yeP, |y-yo|<^ 

for any A G C^{V,R'^) satisfying curl A = B and for any g G i°°({| • — yo| < ^/2}) one has 

TrL2(p)5/((-zV-A)2)-Tri2(R2)g/((-iV-A)2)| <C£2||g||^(/i£ + r^) . (4.3) 

Here (-iV - A)^ and (-iV - A)^ are the Dirichlet realizations of the corresponding operators in 
L'^{V) andL'^{R\), respectively. 

This proposition, combined with our previous analysis about half-space operators leads to 

Corollary 4.2. Given f G S{R) and Af > there is a constant C > with the following property: 
Let q G C^iR) satisfy ||g'||oo + |k"||oo < M and let B G C^{Rl) satisfy \\B\\oo + ||VB||oo < Af. 
Moreover, let ya edV, he (0, 1] and £ G [1, C^^/t.^^/^] ^^^^ ^^^^ g' (2/0,1) = 0, 

sup \q"{y,)\<Mh'/\ (4.4) 
y6'P,|y-yol<^ 

and 

sup |VB(y)| < jj^f B{y)>M-^. (4.5) 

x£V.\y-yo\<e xeP,|y-yo|<f 

Finally, assume that g eC^{{\- -yol < e/2}) satisfies 

sup |,9(y)|+^"' sup \Vg{y)\<M. (4.6) 

ye-P, |y-yo|<^/2 yeV,\y~-yo\<l/2 
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Then for any A G C2(7',]R2) yjH]^ curl A = B one has 



Tr.g/((-zV - A)^) - ^ £ ^ 3iy)bk{B{y), f)B{y) dy 



g{y)sk{B{y)J)^B{^da{y) 

dV 

< Cf [h^l^l^ + + . (4.7) 

Here (— zV — A)^ is the Dirichlet realization of the corresponding operator in L^{V). 

Remark 4.3. Wc emphasize that assumption (j4.4p is weaker than what we need later in our 
application, where we have a similar bound but with the right hand side replaced by Mh. Imposing 
this stronger assumption, however, does not lead to an improvement of (|4.7I) . The reason is that 
the error in (|4.7p comes from Corollary 13.21 and the error coming from straightening the boundary 
in Proposition 14. II is of lower order. 



Proof of Corollary \4-<^ After a translation, if necessary, we may assume that yp = 0. In particular, 
we have q{0) = q'{0) = 0. According to Proposition 14. II we can replace, up to the accuracy we are 
interested in, the trace Trg/((— ■iV — A)^) by Tr^2(Kp *V — A)^). We shall show that the 

latter trace can be computed using Corollarv 13.21 This will prove the corollary, since 



— ^ / g(x)6fe(B(x),/)B(x)dx = — ^ / g{y)bk{B{y),f)B{y)dy 



k=l fe=l ' 

and 



— ^ / g{xi,0)sk{B{xi,0),f)VB{xi,0) dx, - / g{y)sk{B{y)J)^B{^ da{y) 



The last bound follows from the fact that 1 < \/l + q''^ < 1 + Ch£^ on the support of g. This 
bound follows from (|4.2p in Lemma below. 

We now verify the assumptions of Corollarv l3.2l As observed before, A generates the magnetic 
field B{x) ~ B{xi,X2 + q{xi)). Thus the assumption ||-B||oo + ||Vi?||oo < M, together with 
ll'z'lloo ^ M, implies the corresponding assumption for B. In order to prove the local assumption 
on B and g we need the following two inclusions, 

{x e : |x| < A-^L} C {x e : xj + {x2 + qixi)f < L^} C {x e : |x| < AL} (4.8) 

valid for all L £ (0,£] with A ~ 1 + Ch^^^i, where C only depends on M and where we assume 
that h^/^£ < 1. Accepting (|4.8p for the moment, we deduce that the bounds on B and g on the 
sets {y £ V : \y\ < (} and {y G V : \y\ < ^/2}, respectively, will imply corresponding bounds 
for B and g on the sets {x e : |x| < A-^£} and {x € R^^ : |x| < Ai/2}. If h^/^e < (3C)-i, 
where C is the constant entering in the definition of A, then the ratio of the radii of the two disks 
is /2 < (4/3)^/2 = 6 < 1. The assertion now follows from Corollarv 13. 21 taking Remark [23] into 
account. 

We now turn to the proof of (|4.8|) . Assume first that |x| < cL for some c < 1 and L < £. 
Then, in particular, |a;i| < cL and \q{xi)\ < QcL where Q = (A///2)(l + M/A)h'^^'^e from (g^]). 
This implies 

{X2 + q{xi))^ <xl+ 2x2\q{xi)\ + q{xif < xl + 2Qc^L^ + Q^c^L^ 

and therefore x\ + {x2 + q(x\))^ < (?L^{1 + 2Q + Q^). This implies the first part of the assertion 
with any A>l + Q. Conversely, assume that x\ + {x2 + q{xi)Y < for some L < £. Then, as 
before, |(?(a;i)| < QL. This implies X2 < \x2 + q{xi)\ + \q{xi)\ < L{1 + Q) and therefore 

xl < {X2 + q{xi)f + 2x2\q{xi)\ - q{xif < {x2 + qixi)f + 2L^Q{1 + Q) 
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Hence |xp < + 2Q(1 + Q^)). This implies the second part of the assertion with any A > 
1 + 2Q(1 + Q^). Taking into account that h^^^i < 1, we obtain what we claimed in (j4.8p . This 
completes the proof of Corollary 14. 2 1 □ 

Before beginning with the proof of Proposition 14. 1 1 we state and prove 

Lemma 4.4. Assume that q <E C^{R) satisfies q{0) ~ q'{0) = 0, ||<z"||oo < 

sup \q"{yi)\<Mh'^\ 

for some positive M,h,£. Then \q' {yi)\ < for all yi € M. Moreover, for all yi with \yi\ < £ 

\q'{yi)\ < Mil + M/A)h'^^\yi\ and \qiy,)\ < (A//2)(l + M/4)h'/^yl . (4.9) 

Proof. The bound < simply follows by integrating the global bound ||9"||oo < M. 

Similarly, the second bound in (|4.9p follows by integrating the first one. To prove the first inequality 
in (|4.9p . we observe that by assumption we have g'(0) = and \q"iyi)\ < Mh^/"^ for any yi for 
which there is an X2 > such that y\ + {x2 + q{yi)Y — ■ Easy geometric considerations show 
that such xi exists at least when \y^\ < £{1 + (M/2)/ii/2)-i/2 (indeed, the domain Pn {|y| < £} 
lies between the two parabolae X2 — ±{AI/2)h^^'^yf.) Thus for all such yi we have \q'{yi)\ < 
M/ii/2|j/i|. On the other hand, for ^(1 + (M/2)/ii/2)-i/2 < \yi \ < £ wc use in addition the global 
bound \q"\ < M and find 

\q'iyi)\ < Mh^'^{l + (M/2)/ii/2)-i/2£ + M{\yi\ - £{1 + {M I2)h^'^)-^'^) 
< M\yi\ + 1 - (1 + (A//2)/ii/2)-i/2^ 

The first bound in (j4.9p now follows from the elementary inequality 1 — (1 + x)^^^^ < x/2 for 
x>0. □ 

We now give the 

Proof of Proposition \4-l\ The proof is similar to that of Proposition 13.11 and we only sketch the 
major differences. As in the proof of Corollary 14.21 we may assume that yo = 0, so that q{0) = 
q'(0) = 0. We have 

TrL2(p) gfii-tV - A)2) = Tr^.^,^) ~gf{H) , 
where H was introduced before the proposition. We write 

H = Ho + W 

with Ha = (-zV - A)^ and 

W - -(-z9i - Ai)q'{-id2 - A2) - {-182 - A2)q'{-idi - A,) + {-id2 - A2)(q')^{-ld2 - A2) . 

(Strictly speaking, the notation Hq instead of Hq would be more consistent, but we try to keep 
the notation simple at this point.) Similarly as in the proof of Proposition 13.11 we denote by 77 
the characteristic function of the disk {|y| < £}. Using the Helffer-Sjostrand formula we find as 
before that 

Trg(^f{H)-fiHo))^-lJ^dxdy^{z + l) 

X Tr (rXi/o - zy^Ho + 1)-'W{H - z)-'g + g{Ho + 1)" V(F + l)-\H - z)-^fj^ , 

and we need to show that the operators fi{Ho — z)^^{Ho + l)^^W{H~z)^^g and g{Ho + l)^^W{H+ 
1)^^{H — z)^^fj are trace class with a suitable bound on their trace norms. 
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The factors fj{HQ + and {H + l)~^fj arc trace class with norms as in (|3.11l) . Here we also 
use the fact that the support of fj is contained in {|x| < Ai^} with a constant A depending only on 
M, which follows as in the proof of (j4.8p . (At this point we use the assumption that h^^^l < M.) 

Thus it remains to control the operator norm of W{H — z)^^g and g{HQ + We begin 

with the second operator. Again, we decompose W = Wfj + W{1 — fj). 

We now explain how to bound the operator g{HQ + l)~^Wr). We commute all derivatives in 
W to the left and find 

w' = -(-z9i-ii)Ha2-i2V-H92-i2)(-«9i-ii)g' + Ha2-i2)'('7')'-«H52-i2V'. 

Next, we commute one derivative through {Hq + 1)^^. Given j G {1, 2}, let f G {1, 2} \ {j} and 
compute 

(iJo + ^)~\-idj - A,) - {-id, ~ Aj){Ho + 

= (Ho + [-id, - Aj, {-id,, - i.O^] (i/o + 1)-' 

= -i{-iy{H^ + l)-^ (B{-zdy - A,,) + {-idy - Ay)B) {Ho + I)-' . 

Thus 

{Ho + ly^Wfi = - {-idi - Ai){Ho + l)-\-id2 - A2)q'f} - {-id2 - A2XH0 + iy\-idi - ii)g'ry 
+ {-id2 - A2XH0 + l)-\-td2 - A2){q'ffi - i{Ho + l)-H-^d2 - A2)q"v 

- i{Ho + (B{-id2 - A2) + {-id2 - A2)b) {Ho + l)-\-id2 - A2)q'f, 

+ i{Ho + I)-' {B{-idi - ii) + {-idi - {Ho + l)-\-idi - Ai)q'f, 

- z{Ho + [B{-tdi - ii) + {-idi - {Ho + l)-\-td2 - A2){qTfl 

Now we use the fact that \\{-id.j - Aj)°'{Ho + iy^{-idk - AkY^W < 1 for any a,/3 G {0,1} 
and j,k e {1,2}. Since \B\ < M everywhere and since \q'\ < AI{1 + M/4)h^/^e (see (|49|) ) and 
< Mh}/"^ on the support of 77 we conclude that 

\\{Ho + iy^Wfi\\ < Ch^'^l 

with a constant C depending only on M . 

Finally, we bound the operator g{Ho + 1)~^M^(1 — ff). Our goal is to apply Proposition IB. 7[ 
which will give us an exponentially small error with respect to In order to apply this proposition, 
we have to bound the operator norm of 

in terms of M . To do so, we note that W{\ — fi) is given by the same expression as Wfi, except that 
f] is replaced by 1 — r). Hence {Ho + 1)~^ W^(l — ff) has a similar representation as {Ho + 1)~^ Wq but 
with fj replaced by 1 — 77. The required bound now follows from the bounds on B, q' and q" together 
with Proposition ED which states that \\e^^-'^°'^ {-id, - Aj)"{Ho + l)-^{-idk - Ak)'^e-^<--'">m 
is bounded for any a, /3 G {0, 1}, j, k G {1, 2} and \6\ < Sq. 

This concludes the proof of the boundedness of g{Ho + 1)^^W. The bound on W{H — z)^^g 
is derived similarly and we only sketch the major differences. First, it is more convenient to work 
on V instead of and to consider the operator W{H — z)~^g, where H = (— iV — A)^. Now 

W = -{-idi - A^)q'{-td2 - A2) - {-id2 - A2)q'{-idi - Ai) - {-ld2 - A2){qT{~ld2 - A2) 

(that is, only the sign of the last term changes) and we decompose W = ijW + (1 — rj)W . Next, 
we commute all derivatives to the right and then we commute one derivative through {H — z)~^ . 
The argument then is similar to what we have done before. The z-dependence is controlled by the 
resolvent identity in the bound for rjW and by Proposition IB . 71 in the bound for (1 — rj)W . This 
concludes the proof of Proposition 14.11 □ 
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5 Proof of Theorem 11.21 



Instead of working with tlic operator and the fixed domain i7, we will dilate the domain as 
in (|1.12p obtaining = {x € : /i^/^x € fi}. Then will be unitary equivalent with hH^^ 
where Hh = (— iV — A^,)^ is defined in L^{Q,h) with Dirichlet boundary conditions, the new vector 
potential becomes A/i(x) := /i~^/^A(/i-'^/^x) and the new magnetic field is Bh(x) :— B{h^/^x). 
Thus (|1.2p can be rewritten as: 

Trf{Hh) = h-^(Co{f) + h^/^C\{f)+o{h^/^)) asft->0, (5.1) 

where Co(/) and Ci(/) are the coefficients given in the introduction. 

5.1 Partitions and cut-offs 

Our construction will depend on a parameter £, which will later be chosen as an inverse power of 
h. (Indeed, £ = h~^^^ .) However, at this point we prefer to consider £ as a parameter independent 
of h. We will assume throughout that £ > 1, h < 1 and that h^^'^£ < d, where d is chosen so large 
that the domain {x G : dist{x, dil} < lOOd} can be parametrized in terms of a tangential and a 
normal coordinates. This is explained below in more detail. The existence of such d follows from 
the assumption on dfl. 
We define for t > 0: 

Ei{t) {x € I4 : dist{x, 517^} < U} . (5.2) 

This models a 'thin' compact subset of flh, near the boundary, with a volume of order 0{h~^^^£). 
Because we assumed that is small enough, all points ofEi{t), t < 100, have unique projections 
on dilh- We also note that if ti < t2 then S^(ii) C ^^(^2) and: 

dist{Ee{h),nh\Ei{t2)} > {t2 - ti)£ . (5.3) 

The subset r2;i\ 2^(100) will model the bulk region of fth, which is still 'far-away' from the 
boundary. 

Covering the bulk region. 

We choose three functions < 50 < .9o < 9o < 1 with the following properties: 

supp(.go) C r!,, \ S,(100), .go(x) - 1 on 17^ \ (150), (5.4) 

supp(go) C nh\Ei{50), 5o(x) = 1 on ?V\S7(75), (5.5) 

supp(^o) C nh\Ef{10), .^o(x) = 1 on n^\EAW). (5.6) 

There is a constant C such that for all h e (0, 1] and all ^ G [1, dh^^^^] we have: 

max{||i?'3^o||oo, ll-D'^Solloo, P^SolU} < Crl'^l for all |/3| < 2 . (5.7) 

Moreover: 

go9o = 9o, 9ogo = go, dist{supp(Z3go),supp(go)} > ^- (5.8) 



Covering the boundary region. 

We recall our assumption that dfl is a finite union of disjoint regular, simple and closed curves. 
For simplicity, let us restrict ourselves to the simply connected situation. As explained before, we 
assume that h^^^£ is sufficiently small such that 2^(100) can be parametrized using a tangential 



coordinate a living on the torus 



-S^ and a normal coordinate t in (0, 100^). We want to 



divide this strip in curvilinear rectangles of side size proportional to £. We define: 



Me := [/i-^/2^|arj|j e N, 



:= (s - 1) £ , 1 < s < A/f - 1 , an, := ai . 



(5.9) 



20 



Note that the (curvihnear) distance between (7Mi-i and (JMi,{~ ci) is at most 2£ and at least 
£. All other curvilinear distances between consecutive points are £. Note also that due to the 
regularity condition on 90, the Euclidean distances between these points on dflh are of order i. 
Define also: 

C,(ti,f2):={x=(-7,r)GS,(i2): \a~a,\<he}. (5.10) 
We can now find Mi smooth functions {gs} with the following properties: 

supp(5«) C C,(10,150); (5.11) 

Mf-l 

0<gs<l, ^ 5,(x) = lonS,(100); (5.12) 

P'^Ssljoo < C^~"^' foraU|/3|<2, 1 < s < Me - 1 . (5.13) 

In order to use a more compact notation, we denote the set of centers of all curvilinear squares of 
the type Cs{ti,t2) with i?bdy Thus (|5.12p reads as: 

0<5^<1, ^ g^(x) = lonS,(100). (5.14) 

7e-Ebdy 

Another important property is that a point x e (100) can only belong to the support of a finite 
number of g-y^s, and this number is bounded independently of h and £. 

For every 7 € E'bdy we can construct two other smooth functions < g-y < g-y < 1 with the 
properties: 

9i9i=9i, 9i9i=9i-. max{||i:»''.g^||oo, ||£'''.97l|oo, ||£'''.g7l|oo} < C£~l'^l for aU |/3| < 2, 

(5.15) 

dist{supp(Dg7), supp(57)} > £ , diam{supp(5-y)} < . (5.16) 



Patching the two regions together. 

It is clear from the construction of go and the boundary covering that they can be made compatible 
with the supplementary condition that they generate a partition of unity of the whole dilated 
domain fJ^. Thus we can impose: 

5o(x)+ Y. •97W = 1, Vxe^!,,. (5.17) 

76-Ebdy 

5.2 Local model operators 

Let us define 

Ar^x) :=.^o(x)A,(x). (5.18) 

We can define an operator Hf^"^^^ := (— iV — A^"''^)^ defined in _L^(M-^) and see that due to the 
support conditions we have the identity: 

goA}:""^ = , iJ^^'^go - i?/.5o . (5.19) 
Moreover, for every 7 G -Bbdy we define: 

A|7)(x) :=^^(x)A,(x). (5.20) 

For all sufficiently small h^^^(, the support of each g^ contains a piece of 9ri/i having a length of 
order 0{£) and is, after a translation and rotation, included in a perturbed half-plane P^'^'^ of the 
form (|4.ip . Here q^'^'> coincides with the local parametrization of dflh- We now define the operator 
Hj^'' ~ (— iV — ^'h^)^ with Dirichlet boundary conditions in L'^iV^''^). Again due to support 
considerations, 

-g^A'l^^ ^-g^An. h[^^ ~gy ^ Hh~g^ . (5.21) 
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5.3 An approximation for f{Hh) 

The following lemma tells us that, up to a controlled error, we can compute Tr f{Hh) by computing 
Tr f{Hf^''^^)go and Tr f{H\^^)g^ for 7 S iJbdy. We recall that d is the upper bound on h^/'^i, which 
is needed in order to introduce boundary coordinates. 

Lemma 5.1. Given f € S{M.) and N > 0, there is a constant C > such that for all h e (0, 1] 
and ^ € [1, dh~^^'^] one has 



7e-Ebdy 



(5.22) 



Presumably, the order h^'^/^ in (|5.22p is not optimal, but it is sufficient for our purposes. 



Proof. Step 1. An approximate resolvent 

If z e C with Imz 7^ 0, we define the operator: 

7e-Ebdy 

One can prove that the range of Sh{z) is in the domain of Hh and we have: 

{Hh-z)Shiz) = l + Vniz), 

with 

Vn{z) := {-2*(Vgo) • (-*V - A^^'^) _ (a^^)} (i/buik „ ^yi^^ 

+ J2 {-2*(V,9^).(-zV-A|7))-(A.9,)}(ff(^)-z)-\g,. 

76-Ebdy 



(5.23) 



(5.24) 



(5.25) 



In order to obtain (|5.24p we observe that Hh can be replaced by the local Hamiltonians on the 
support of the corresponding cut-off functions, see (|5.19p and (|5.2ip . The operator Vh contains 
only terms coming from commuting the local Hamiltonians with cut-off functions. The identity 
operator on the right hand side of (|5.24p appears after the use of (|5.17p . 
It follows from (|5.24p that 

(Hh - z)-' = Sh{z) - [Hh - z)-^Vh{z). (5.26) 

Standard estimates (see (jB.Sp ) yield for any h G (0, 1] the existence of a constant C{h), such that 

lRezp/2 

\mz)\\<C{h)L—^, WzeV. 
|ini z| 

The Helffcr-Sjostrand formula p.6p shows that for any almost analytic extension fa,N of / satis- 
fying ([2|) with > 3 one has: 

/(ff/0 = 5o/(i^r')5o + E hf{Ht'^)9-i-- I ^k^{Hu~z)-^V,,{z)dxdy. (5.27) 



Hence (|5.22p will follow if we can show that for any A^i > there is an 7V2 > and a C > such 
that for any h € (0, 1] and ^ € [1, one has 



df N 

— ^^{Hh - zy^Vh{z) dxdy 



(5.28) 



Bi{L^{Qh)) 
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Step 2. Preliminary estimates 
Wc show that {Hh — z)~^ is a Hilbcrt-Schmidt operator. First, for z = — 1 we use the diamagnetic 
incquahty in order to bound the absolute value of the integral kernel of {Hh + with the 
(positive) integral kernel of (—A + 1)^^. Thus we obtain: 



\iHh + l) 



IS2{L2(Oh)) 



(5.29) 



For general z we use the resolvent equation 

{Hh - z)-^ = {Hh + + {z + l){Hh - zr^Hh + 
in order to obtain the existence of C > such that 

(Rez) 



(5.30) 



for all h and z £ D. 

Next, we consider the operators Vh{z) and prove that there arc constants C and 6 such that 
for all z e all ft. G (0, 1] and all £ e [1, d/i^^/^j ^isls 



\\Vh{z)\\ < C/i-3/4£-i/2(I^exp{-5^|Imz|/(Rez)} 



and 



\\Vh{z)-Vh{-l)\\s,^L.^n,)-, 



llmzl 



Imz 



■exp{-d£\lmz\/{Rez)}. 



(5.31) 



(5.32) 



In the following z will always be restricted to T>. In order to prove (|5.3ip . let us take a generic 
term of Vh{z), that is (Vg-y) • (— iV — aI^^){hI^^ — z)~^g'y. We can cover the supports of Vg^ 
and g-y with approximately i'^ disjoint squares of diameter 1. Denote one such square covering the 
support of Vgj by Ai, and one such square covering the support of g-y by A2. Let Xj: J = 1)2 
denote their characteristic functions. Due to the condition (|5.16p . the distance between Ai and 
A2 is at least of order £. 

Now nj^"^ plays the role of Hjj in (|B.13p . and (Vg-y) • (— iV — ASf^"^) is W in the same estimate. 
A consequence of Proposition IB . 71 and of (|5.15p is the estimate: 



Xi(V.g^) • (-*V - a|7))(f(^) - z)-'x2 



< Cr^^^exp{-Se\lmz\/{Rez)}, 
Im z 



where the constants C and S are independent of h, 7, z and the positions of xi and X2- 

A similar estimate holds for xi{^9-y){Hl^'^ — z)~^X2- Now let Xi'^\ fc = 1, • • • , A', be a family 
of characteristic functions corresponding to disjoint squares of diameter one covering supp(Vgy), 
and similarly X2 \ ^ = Ij ■ • ■ 7 covering supp^.^. Then for all 5*, $ € L'^{flh) 

{-2^{V~g,) ■ {-^V a|7') - (A.g,)} (4^' - z)-\g, $) 
< E {-2^(V.g,) • (-zV - A^^) - (Aff,)}(4''' - zr'g,}x^i^^) 

I (Re z 



k,l 

< ce 



|Im z 



■exp{-J€|Imz|/(Rez)}E Xi * X2 * 



k,l 



Since K and L are bounded by a constant times we have J2k 
inequality and similarly for From this we conclude that 



< C£||*|| by the Schwarz 



z) ^g-f 



£jRez) 
llmzl 



'|Imz|/(Rez) 
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Summing over 7 and using the faet that number of 7's which arc non-zero at a given point 
is uniformly bounded in h and £ we arrive at (j5.3ip . The estimate for the term containing 
^buik jg completely analogous, with the difference that the number of unit squares needed to 
cover the support of go is of order whereas the number of unit squares needed to cover the 
support of Vgo is of order This explains the factor h~^^^£~^/'^ in (|5.3ip . The factor 

£ does not appear in (|5.31|) because of our assumption that h^/^£ < d. (It is the bound on 
{-2i(Vgo) ■ (-iV - Al""^^) - {Ago)} {H^"-'^^ - zy^go, which is not optimal. One should take into 
account the decay between the supports of xi''^ ^nd similarly as we did in the proof of 
Proposition 13. in 

In order to prove (j5.32p we proceed similarly as before by controlling the Hilbert-Schmidt norm 
of the operator: 

Xi(V.g,) • (-*V - A^:^){{hI:^^ - z)-' - {H^^^ + 1)-']X2 

= {z + l)xi(V.g,) • (-*V - a(-^')(4^) + - z)-^X2 

with xi and X2 as above. Now using the result (|B.18|) of Proposition IB . 71 we arrive at: 



< (^^-li5£^g-«|Imz|/(Roz) 

Bo ~ \lmz\ 



where the constants C and 5 are independent of /i, 7, z and the positions of xi and X2- Hence, 
with the same notation as above. 



= E ||xf '(V5.) • (-^V - A^:^){Hi^^ + irHH^:^ zy^g^xf 

k,l 

< Ce ^^""''Z exp{-2J£|Img|/(Rez)} 
|Im 

Here we used that both the k and the / sum contain 0{£'^) terms. Summing over 7 and using 
again that the overlap of the covering is uniformly bounded, we arrive at an estimate as in ()5.32p . 
The term with II]^^^ can be treated in a similar way, where we note that we need about h~^^ unit 
squares to cover the support of go and h~^/'^£ to cover that of Vgo- (Again, the argument can be 
improved at this point.) 
Step 3. Proof of ({5^281) . 

Due to the presence of the momentum operator in the expression of V/i(z), we cannot hope to 
prove that Vh{z) is Hilbert-Schmidt. We resolve this problem in a similar way as we did in the 
proof of Proposition 13. 11 namely by replacing the operator (H^ — z)~^Vh{z) under the integral in 
([Qg]) by the difference 

(Hh - z)-'Vhiz) - {Hh + ir'Vhi-l) 
without changing the integral. In this way, we obtain a sum of two terms: 

[{Hh - - {Hh + ly^] Vh{z) + {Hh + I)-' [Vh{z) - Vh{-1)] . (5.33) 

In the first term we can apply the resolvent formula and arrive at the operator: 

{z + l){Hh - z)-\Hh + ir'Vh{z). (5.34) 

By (|5.29p and (|5.30p the product {Hh — z)^^{Hh + 1)^^ is a trace class operator with trace norm 
bounded by C h"-^ (Re z) /\lu].z\. The bound ()5.3ip on the operator norm of Vh{z) gives 



\{z + l){Hh - z)-\Hh + 1)- V„(z)|L < Ch-'IH-^/^^^^ exp{-<5^|Imz|/(Rez)} 
I wax Imzr 
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which, in view of the inequahty e ^ < (N/e)^ , can be further bounded as foUows: 

\\iz + l)iH, - z)-\H, + < C^fe"V4^^A^-i/2 (R^y^ _ (5.35) 

Thus we obtained a power-hke decay in h at the price of a more singular behavior in Imz. But 
choosing N large enough in (j2.5p allows us to kill the singularity in Im z and - at the same time 
- to turn the growth in Re z into something integrable. Therefore we are done with the first term 
of (1533)) . 

The second term in (|5.33p is bounded using (|5.29p and (|5.32p : 

II (if, + [Vhiz) - Vh{-l)]\\g < Ch-^^H-'/^^-^^ exp{-<5£|Imz|/(Rez)} . 
II iiai |imz| 

As before, this implies 

WiH, + D- [V,iz) - V^.(-l)]||,^ < C.h-^^'r--^^'^^^ 
which completes the proof. □ 

5.4 Completing the proof of Theorem 11.21 

The bound from Lemma [01 reduces the study of Tr{/(_ff/i)} to that of 

TbuikW :=Tr{/(i/rV} 

and 

76-Ebdy 

Indeed, Lemma [01 leads to the estimate: 

Tr{/(H,0} = TbuikW + Tbdy(M + (5.36) 
for any > 1. Let us briefly mention that by writing 

we can apply Proposition [B3] and conclude that f{H^i^^)g^ is trace class. The same is true for the 
bulk term. 

The bulk contribution. 

Since the operator /(-ffj;'""^)5o is trace class and has a jointly continuous integral kernel (see 
Appendix ICj). we can compute its trace as Ti[h) = f{H}^^^^){^,^)go{^)dx. The bulk operator 
^buik constructed in such a way that on the support of go the vector potential A)j""^ generates 
the magnetic field Bh{x), which obeys p.ip . Now putting together p.2p and p.3p we obtain: 

00 „ 

Tbuik(/i) = (27r)-i ^ / bk{Bh{y),f)Bh{y)go{y)dy + ON{\snppgo\{h + r^)) 

k=i •'^t- 
00 „ 

= (27r/i)'i V / 6fc(B(x), /)B(x)5o(/i-'/'x) dx + On{1 + h-'r'') . (5.37) 



The boundary contribution. 

The operators Hj^^ are constructed in such a way that the assumptions of Corollary 14.21 are 
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satisfied, with a constant AI independent of 7. The role of g will be played by g^. On the support 
of g^ we have introduce a local parametrization g^'*'^ of fi/j. Since Vlh is a dilation of f2 by a 
factor the functions /ii/2|^(7)'| and are bounded on the support of g^, uniformly 

in h, £ and 7. This is, of course, stronger than (14.41) . Thus there is a C such that if ft, e (0, 1], 
i e [1, C"^/i~^/^] and 7 e E'bdy one has 

00 „ 

Tr{/(if(^))5,}--n-E / 57(/^~'/'x)6,(i?(x),/)i?(x)dx 
fe=i 

-TTTTlE / 57(/^"''"x)sfc(i?(x),/)v/BWda(x) 

CO „ 

fe=i "'^^h 

Since there are about h^^^^£~^ elements in E'bdy we obtain 
00 „ 

Tbdy(M-^E / 9yih-'/^^)b,iBi^)J)B{^)d^ 

^ k=l •'^ 

— Y. / 57(/*"'/'x)sfc(i?(x),/)/BWda( 



27r/ii 



x) 



< Ch-^''^£{h}''^£^ + M'^ + r ^) . 

(5.38) 



In order to optimize the error, we now choose £ = h^^^^. (Note that with this choice h^^^£ is 
arbitrarily small if h is small.) With this choice the right side of (|5.38p is bounded by a constant 
times 

Finally, by choosing N in (|5.36p and (|5.37p sufficiently large and recalling ()5.17p . we arrive at 
Tr f{Hh) = h-' (Co(/) + h^'^CiU) + 0(ft5/8)) as /w . 
This completes the proof of Theorem 11.21 



6 Proof of Corollary 11.31 

This corollary is a rather direct consequence of Theorem ll.2l We will work with the dilated domain, 
see (|1.12p . Clearly, because Lh is unitarily equivalent with hHh we obtain that N{hE,Lh) will 
equal the number of eigenvalues of (including multiplicities) which are less or equal than E. 
In order to simplify things we assume that K = 1, which means that B-^g^ < E < 3i?min- 

Theorem 11.21 can now be applied for a function / £ (M) which equals 1 on the inter- 
val [Bmin — e, Sniax + e] (hcrc e > is chosen small enough) and whose support is included in 
(— oo, 3i?min ~ £]• In this case only the first Landau level will contribute to the bulk term, and we 
can thus state the following straightforward consequence of our main theorem: 



lim 

h\0 



TT{f{H„)}-— y^^S(x)dx 



^ I da(x) r dt^w^) 

Jon Ja 



f{B{^))e,{0)\Mt,0\ dC-l 
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Now denote by /_ any function as above whose support is included in {—oo,E], and by /_|- any 
function as above which equals 1 on {^oo,E]. By a straightforward monotonicity argument we 
have: 



1 

2^ 



da(x) / dt^B{x)\ / /_(S(x)ei(0)|^i(i,Or'^e-l 



< liminf /i^/^ 



< limsup/i^/^ 



Nh{hE)-^ f B(x)dx 



< 



^ [ da(x) r dt^/m^i [ /+(B(x)ei(0)|V'i(i,Ol'rfe-l| • (6.2) 

At this moment we can take the supremum over /„ and the infimum over /_|_ and show that they are 
equal. The main idea is to show that the Lebesgue measure of the set of ^'s near i3(x)ei(^) = £' is 
small, uniformly in x G dfl. In fact, using Hadamard's identity (|A.9p one can prove that uniformly 
in X g dH.: 

lim\{^eM: E -S < B(x)ei(^) < E + S}\ = . 

»5\,0 

The last ingredient is to use that ipi has norm one. The proof is over. □ 



A Analysis of the one-dimensional model operator 

A.l Eigenvalues and eigenfunctions of the model operator 

Recall that the eigenfunctions (f)k and ■(/'^(cot,^) were introduced in the introduction. Here we 
collect some information on these functions as well as on the corresponding eigenvalues. 

Lemma A.l. There are constants a > and C > such that, for all k > 1, 



\Mt,0\<C^/e,.{0 + l, \Mt)\< 


C Vk, t>0,^eR, 


(A.1) 


efe(C) > (2fc-l), VeeM, 




(A.2) 






(A.3) 


\Mt)\<Ce-"'\ \t\>Ck^, 




(A.4) 


|efc(e)-(2A:~l)| <Ce-"l«^ 


-Cki, 


(A.5) 




t > and ^ < -Cki . 


(A.6) 



Remark A.2. The main observation is that tpk lives in the classical region + < efe(^). Hence 
outside a sufficiently large neighborhood of this classical region, the eigenfunction should decay 
exponentially. The second point is that if (the boundary point) is far from the classical region, 
then (/)fe(i + 1) furnishes an excellent approximation of ipk in large zones containing the classical 
region but avoiding a neighborhood of 0. This is what we now have to control by using explicit 
Agmon estimates. 

Proof of Lemma \A.1[ Let us first prove (|A.ip . We recall that V'fc and (f>k have norm equal to 
1. One immediately obtains after an integration by parts that 

Using the Sobolev embedding, we get the existence of C > such that, for any k > 1 and ^ G K. : 

snp\Mt,0\ < C^l + ekiO . (A.7) 

Similarly, we get 

snplMs)] < CVk . (A.8) 

s6R 



27 



The two inequahties (|A.2p and (jA.3[) are easy by monotonicity of the Dirichlet problem with 
respect to the domain. Let us also note for the record the Hadamard formula applied to our 
Dirichlet problem: 

^e,(0=^UO,ef . (A.9) 
Let us show that (|A.4[) is a consequence of an Agmon estimate using a weight e**-"^ with 

- /3y , < /3 < 1 . 

Since {D^ + — {2k — l)}(j)k{s) = , we obtain the identity: 
or 

We^Mhm = Jj2k - (1 - l3')s^}e^''^^Ul{s)ds < 2ke^ , (A.IO) 

where the inequality was obtained by taking a rough upper bound of the factor multiplying 4>k{s) 
on the region where this factor is nonnegative. Hence the Sobolev embedding provides a constant 
C > such that 

\Ms)\ < Cke^ e-'^"'/^ s e M. (A.ll) 
Now, for < /3 < 1, there exists such that 

Cfce^e"'^"'/^ < 1, \s\>CpVk, 

and (jA.4l) follows by choosing for example a = /3/4. 

Also note that ()A.4p implies the existence of some smaller a and possibly larger C such that 

\Mt + 0\<e-<^'^^^"+^") ,yts.t. |t + ei>y>Cy^. (A.12) 
Let us show (|A.5[) . Due to the lower bound (|A.2p . we only need to show the upper bound 

efc(^) < (2fc- 1) + Ce-"I«I' , fore<-CVfc. (A.13) 
For a given a, we can always find C such that Ce~"'"^ < i. In particular, we have 

(2j- - 1) + Ce-"l«l' < (2j + 1) - i , ^<~C, j>l. 

We will show that, for every 1 < j < k and for every ^ < —CVk, there is in an interval of half 
length Ce~"^ centered at (2j — 1) at least one eigenvalue ee{^) of the Dirichlet realization of 
-Dj + (i + C)^ in R+. Since (jA.2p implies that we cannot have more than k eigenvalues in the 
energy interval ]1, {2k + each such constructed ei{^) must be alone in its own interval. 

We take x{'t)4'j{t + the j-th quasimode, where < x < 1 is a C°° cut-off function equal 
to 1 on [1, +oo) and vanishing in (— oo, 0]. The control (jA.4[) (which also implies the control) 
yields 

||(i?^ + (t + e)')x</',(-+e)|L.(R^) <C^e-<' and ||x0,(- + C)|U^(r^) > 1 - Ce'^' (A.14) 

for ^ < —C\fk. The spectral theorem now leads to the conclusion. Hence, we have shown that 
there exist constants C and a > (after possibly increasing the old C and lowering the old a), 
such that 

|ej - (2j - 1)1 < Ce-"«' , < -C^fk, 1 < j < k. (A.15) 
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The bounds (jA.14[) and the spectral theorem also imply that 

Uj-XV3{-+0\\LHM^)<Ce-"i\y^<-CVk, l<:i<k. (A.16) 

We now turn to the proof of (|A.6p . Once we know that Ck ^ 2k ~ 1 under the condition 
^ < —C\/k, what we have done for the Gaussian decay of ipk can be also done for any real 
eigenfunction ipk corresponding to ek- Thus for some (3 > we have 

\Mt)\ < C'/3e-'^l*+«l' , Vt > s.t. \t + ^\> CpVk, C < -CVk. (A.17) 

We note that under the condition ^ < —C\fk, this implies also 

\i>kit)\ < C^e-^(l*+«l'+«'' ,yt>0 s.t. |t + > -| > CVk/2 . (A.18) 
The bound (|A.4p together with (|A.16j) implies that 

Uk-M-+mLHR^)<Ce-"^\ ii^<-CVk. (A.19) 
To have a pointwise estimate, we integrate by parts and observe the following inequality: 

yyk~'p'k{-+of + \\it+oiA-M-+of 

< {2k - 1)11,^, -ipk{-+ Of + '2{ek - {2k - 1)) + |V4(0) - v'k{i)\ ■ (A.20) 

The right hand side is exponentially small like exp(— a^^) if we control ip'k{0 and V'fc(O). This will 
come from a control of both of them in iJ^(]R+), through {D^ + s^) with the appropriate boundary 
conditions. Hence we have 

\\Vk\\H^(R) + ||V^fe||ff2(R_^) < C{{2k - 1) + ek} . 

Thus 

sup \Mt) -^k{t + 0\< Ce-"«' , if e < -CVk . (A.21) 
t>o 



We note that the last inequality implies: 



\Mt)-^k{t + 0\<Ce-^'-^"+^'+^^"K if^<-C\/fcand|i + e|<|i. (A.22) 

We are now able to complete the proof of (jA.6P . This is simply a consequence of (jA.12p , (jA.lSP 
(when < t < ^ and t > ) and of (|A.22[) when ^ <t < . More precisely, we have shown 
the existence of C and a > such that 

iVfe W - Mt + 01 < Ce~^(«'+(*+«)') , if ^ < -CVk and i > . (A.23) 

This finishes the proof of Lemma lA.ll □ 

A. 2 Proof of Lemma 11.11 

Replacing f{E) by f{B^^E), we may assume in the proof that B = 1. Since (pk is normalized we 
have 

' \cPk{t + £,)?di = l, 



and therefore we can rewrite sa:(1, /) in the more symmetric form 

Sk{lJ)^ (^Jjf{ek{0)\Mt:0\'- f{'2k-l)\Mt + 0\') rfc) dt. (A.24) 
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In this form, things become elearer. We can use the estimates from Lemma lA.ll on the behav- 
ior of the eigenvalues efe(^) and eigenfunctions ipkit,^) to show that Sfc(l,/) is weU-defined and 
summable with respect to k. Before we do this, we note that once the convergence of the inte- 
grals in (|A.24[) is proved we can apply Fubini's theorem and interchange the order of integration. 
Since ipk is normalized, this yields the equivalent formula p.lO[) for /). Also the continuous 
differentiability with respect to B > is a consequence of our bounds. 

We now prove that the integral ()A.24|) converges and is summable with respect to k. We split 
the integral with respect to ^ in four regions: (I) ^ > CVk, (II) < e < CVk, (III) -CVk > e < 0, 
and (IV) ^ < —Cy/k, where C is the constant provided by Lemma [A. II 

• We start with region (I). We use the fact that / is rapidly decaying and get that, for any 
A^, there is a Cn > such that 

\fm<CNii + E)-^^ . 

Then we use (jA.Sp and (jA.2[) and write: 

\f[ekm<CNk-^{i+e)-'' ■ 

(Cjv here denotes a generic A'^-dependent constant, which may vary from line to line.) This 
when multiplied by \fpk{t,0\'^ is absolutely integrable in ^ and t, and the resulting integral 
is summable in k ii N > 1. 
For the second term we use (jA.4|) : 

OO /'OO /"OO />oo 

dt / dei/(2fc - 1)1 \Mt + 01' < C|/(2fc - 1)1 / dt e-"*%-"«' < C^fc-^ , 

JcVk Jo Jo 

which again is summable if A^ > 1. 

• We continue with region (II). The first term, /(efc(^))|'(/'fc(i, treated as in region (I). 
For the second term we use the fact that 1/)^ is normalized to get 

dt / dei/(2fc - 1)1 \Mt + of < C\f{2k - l)\Vk < CNk-^+^^^ 
Jo 

which is summable for A^ > 3/2. 

In region (III) we use (|A.2[) and see that both f{ek{0) ^-nd f{2k — 1) are bounded from 
above by a constant times k~^ . The integral with respect to ^ gives just an extra y/k, while 
the integral in t is taken care of by Tpk and 1/)^ . 

Finally, in region (III) we split the integrand in the following way: 

fiek{0)\Mt,0\''fi'2k-l)\Mt + 0\' 

= ifiekiO) - f{2k - 1)) \Mt, + /(2fc - 1) (Mt: + Mt + 0) {Mt, - Mt + 0) ■ 

We treat the two terms on the right side separately. 

Let us start with (/(efc(^)) — f{2k — 1)) \ipk{t,0\'^ . Since we are in the region where (jA.Sp 
applies, we can, by using the mean value theorem and the fast decay of /', bound the double 
integral by: 

C^fc-^ / dt / d^e-"«' \Mt,Of < CNk-"" . 



This takes care of the first term on the right side. For the second term we use (|A.1|) and 
obtain: 

i/(2fc - 1) (^fc(i,c) +Mt+ m < CNk-^^+^f' . 

Then (jA.6|) provides the absolute integrability of i^k{t,£_) — (f>k{t + ^ on the whole region, 
whose integral is uniformly bounded in k. 
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This concludes the proof of Lemma 11.11 

Lemma A. 3. Let f he a Schwartz function on [0, oo). For B > the sum 



k=l 



{f{Be,{0)\Mt,0\^-f(B{2k-l))) 



dt 



converges. 



Proof. The proof is similar to that of Lemma ll.ll Again we can multiply f{B{2k—l)) by 

without changing the integral. We split the ^-integral into the same four regions. The only new 

ingredient is the bound 

oo j'OO POO 

t\Mt,0\^dt= (t + ^)\Mt,0\'dt-d \Mt,0\'dt 

Jo Jo 

•oo \ 1/2 

< I {t + 0^\Mt,0\^dtj +max{-C,0} 



< Ve/c(0 +max{-C,0} 



and the corresponding bound 



t\Mt + 0\^dt< V2fc-l + max{-^,0}. 



We omit the details. 



□ 



B Uniform exponential estimates 
B.l Integral bounds on the Green's function 

In this subsection we assume that A C is an open set and that A = {Ai, A2) G L-^'^°^{A, M^) is 
a magnetic vector potential. We denote by Hd := (— jV — A)^ the operator corresponding to the 
closure of the quadratic form |(— *V — A)ip\'^ dx initially defined on C^(A). It is remarkable 
that the results in this subsection do not require any smoothness of A or OA. We emphasize the 
crucial fact that all constants below can be chosen to be independent of A and A. 
Throughout this section we abbreviate 



V:^{zeC: < |Imz| < 1} 



and 



r := \/ (Re z) , 77 |Imz| 

for z G C. Moreover, if xq is a fixed point in A and a G R, we denote by e"^'~^°'^ the multiplication 

operator with the function x h- > e" v'^^^'^^. The following proposition is the key to all the other 
estimates in the section. 

Proposition B.l. There is a Sq > such that for all xq G A, all z G D and all 5 G (0,(5o] one 
has 



< 



(B.l) 



< 6— , 



(B.2) 
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Its proof relics on the following elementary 
Lemma B.2. For any 2 G C and any j e {1, 2} we have 



||(-iVj - Aj){HD - z)-^\\ < Vl/|Iniz| +max{Rez,0}/|Imz|2 . (B.3) 

Indeed, this is an immediate consequence of the following identity, valid for every ip £ L'^{A), 
2 

J2 - Aj){HD - z)-^^pf = Rc{{Hd - +Rez \\{Hd - ■ 

Proof of Proposition \B.1\ Since multiplication with the exponential weight is an unbounded op- 
erator and does not leave invariant the domain of Hd, we must work with regularized weights of 
the type e"'''^'') where: 

/^.(x):=^^^, .>0. (B.4) 
V^e|x-xoP + 1 

One should note that is bounded (though not uniformly in e), while derivatives of any order 
of /le are not only bounded, but uniformly bounded in e and Xq. In particular, e"''^^''^ leaves 
invariant the domain of Hq. 

For any s > the Combes-Thomas rotation [B] reads 

e/^'{HD - z)e-^''' = Hd - z + 2isVh, ■ [-iV - A] + sAh, - s^lV/i.p. 

Since first and second derivatives of /i^ are uniformly bounded in e and xq, the bound ()B.3P yields 
for s = 5rj/r with < 5 < 1 the estimate 

||{2isV/ie • [-iV - A] + sAh, - s2|V/i,p} {Hd - z)-^\\ < CS . 

Here C is a constant independent of z, S, e and xq. Thus if we choose 6 < 1/(2C) wc may write 

e''''{HD- zy^e-"'"' = {Hd-z)-^{1+ [2isWh, ■ [-iW - A] + sAh, - s^\Wh,\^] {Hd - z)-^y^ 

(B.5) 

and deduce \\e'''''{HD - z)-^e-'''^\\ < 2/t]. For Vi e icomp(A) and i'2 G ^^(A) with = 
1 1 '02 1 1 = 1 this implies 

\{e^^-^"Uji,{HD ~ z)-'e-^^-^"^4'2)\ = lim|(V'i,e^''=(i?i5-^)"'e"^"^V'2)| < 

«N0 77 

By a density argument we see that {Hd — z)~^e~~^^'^^'>'>Tp2 is in the domain of e^^'^''"^ and 
then (jB.ip follows. 

Regarding (|B.2I) . one uses {Hd + 1){Hd - z)~^ = I + {z + 1){Hd - z)"^ on the right hand 
side of (|B.5I) . then after a short density argument one employs (jB.ip and the proof is over. □ 



We also need the following small variation of Proposition lB.il 

Proposition B.3. There are constants 6q > and C > such that for all xq G A, all j, fc G {1, 2}, 
all a,/3G {0, 1} and all d G (0, Sq] one has 

Proof. We use the same notation as in the previous proof. Since 

e'''''{-iWj - Aj)°' = {-iWj - Aj + isVjhXe'''' 
and {-iVk - Akfe-'''^ = e"^''^ {-iVk - Ak + isVtKf , 
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where V/ig is uniformly bounded, wc only need to show that for all sufHciently small s, the operator 
(ffD + l)^^^e*^=(£f_D + l)"^e"^''«(i7D + l)^^^ is bounded uniformly in e e (0, 1]. As in the previous 
proof we can shown that er"'^' {Hd + 1)6"'^' ^ Hd + I + W with \\W{Hd + l^^'^W < C{\s\ + s^) 
uniformly in e £ (0, 1]. Thus for \s\ < min{(4C)"^, (4C)"^/^} wc have 

and 11(1 + {Hd + 1)-^'^W{Hd + < 2. □ 



In the following we denote by Q{-,-,z) the integral kernel of the operator {Hd — z) ^ for 
z G C\ [0,oo). 

Corollary B.4. There are constants 5 > Q and C > such that for all z Cz V and for a/Z x G A 

/ dx' e^<^-^'> |g(x, x', < c4 ■ (B.7) 
Ja T 



Proof. Since CJ(x, x',z) = C/(x',x, z) it is equivalent to prove that there are constants 5 > and 
C > such that for all z E V and for all x' G A 

/ dxe^<''-'''>|g(x,x',z)p <c4- (B.8) 
Ja T 

Due to the diamagnetic inequality, |tj(x, x', — 1)| is pointwise bounded by the kernel of (— A^i + 
1)""'^, which in its turn is bounded by the kernel of the resolvent of the Laplace operator defined 
in the whole of . From the properties of the latter kernel we deduce that 

sup/" rfx |g(x,x', -1)1^ oo. (B.9) 

Moreover, by the resolvent identity 

g(x,x',z) = g(x,x',-l) + (z + l) / dyg(x,y,z)g(y,x',-l). (B.IO) 

Hence 



/ dxe^<^-^'>|g(x,x',z)p<2 / dxe^<^-^'>|g(x,x',-l)|2 

J A J A 

+ 2|z + l|2 / dxe^^^-'^') / dy5(x,y,z)5(y,x',-l) 

J A J A 

By (|B.9|) the first integral on the right side is finite if 5 < 1/2 (because then also 26ri/r < 1). If, in 
addition, we choose S not larger than the Sq in Proposition lB.il then by (|B.1[) the second integral 
on the right side is bounded by 

2^^ 



) |^dxe^<— ')|tj(x,x',-l)|\ 



.V/ Ja 

which again by (|B.9|) is finite. This proves (jB.Sp . □ 

As an application of Corollarv lB.4l wc derive a useful trace class criterion. 

Proposition B.5. There is a constant C > such that for any set Ai C A o/ diameter < 1 the 
operator xi{Hd — *)~^; Xi being the characteristic function of Ki, is trace class and satisfies 
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Proof. Wc fix xq G Ai and write: 

(B.ll) 

The operator ^e^^^"^'^ is bounded by since Ai has diameter one. The last two operators are 
both Hilbert-Schniidt because their integral kernels are square integrable by Corollary IB. 4| and 
their Hilbert-Schmidt norms can be bounded independently of xq. □ 

We now turn to the proper topic of this subsection, namely exponential decay estimates. By 
this we mean bounds in various norms on operators of the form xiW{Hd ~ z)~^X27 where xi and 
X2 are characteristic functions of disjoint sets. 

Proposition B.6. There is a constant C > such that for any disjoint sets Ai, A2 C A of 
diameter < 1 and of distance d := dist(Ai,A2) > 1 and for any z € T> the operator xi{Hd~z)~^X2, 
Xj being the characteristic function of Aj , is Hilbert-Schmidt and satisfies 

2 

WXliHo - zr'x2\\B,(L^A)) < 

Proof. Fix xo G Ai. We can write: 

X^iHn z)-\, = xie-l--l {e-^l--l(F,, - z)-' e'^^-^-^ x^} e--l---lx2. (B.12) 



The operator xie~i^'''°~'' is bounded by e"?, while e~~^'''''~ 'x2 is bounded by e~~^'^. The operator 
in the middle is Hilbert-Schmidt because its integral kernel is square integrable, as can be easily 
inferred from (jB.7p . □ 

Now let us consider the perturbed case. Assume that W is an operator relatively bounded to 
Ho with a relative bound less than one, which obeys the following condition: 



. An I 



cs 



sup sup 

xoGA \a\<S 



_ ,■^-lp-"(•-xo) 



W{HD-iye 



< 00. 



(B.13) 



The example we have in mind is = —idj — aj; a quick commutation shows that this operator 
has the property 



sup sup 

xqGA |a|<5 



< 00, 



but this is not quite as in (jB.lSP : in order to show that (|B.13P holds, we write 

ga(-xo>yy(^^ - i)-^e-°'^-''"^ = (B.14) 
|e°<-^") M^e-"<-^"> {Hd - ^)"^ } { (Hd - i)e"<-'"'^ {Ho - i)-^e-"<-'">^ } 



and now we can invoke (|B.1|) and (|B.2|) to see that C5 < 00. 
Returning to general operators W we state: 

Proposition B.7. There are constants C and 5 > such that for all Ai, A2 C A as in Proposi- 
tion \B.(^ and all z G T> one has 



\xiWiHD + l)-'x2\\<csCe-''-''-'^ 
IxiWiHo ~ z 



, r^ Sr,{d-1) 

-'x2\\<csC-er^^., 



\xiW{Hd - i)-\HD - z)-\2\\b, < cs C—e-'-^^. 

■n 



\xiW{Hd + 1)-\Hd - z)-\2\\b, < csC—e ^ 



5T?(d-2) 



(B.15) 
(B.16) 

(B.17) 

(B.18) 



Here cs is the constant from (jB.lSp . 
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Proof. Denote by Ti(z) ;= {Hd - - z)-^ = 1 + (z - i)iHD - z)"^ and write 

XiW{Hd - z)-^X2 = XiW{Hd - i)-iri(z)x2. 

For the proof of (|B.16|) we insert exponentials as in (|B.12p and use (|B.13p and (IB.1[) . For the 
proof of (jB.17p we use the same idea of sphtting the exponential decay as in (jB.llI) . Namely, if 
we choose xq in the support of xi we can write: 

XiW{Hd - i)-\Hd - z)-\2 = xe^l"°-l {e"^!''"-! W^(iJz5 - ^)-le^l''°-l } (B.19) 

.|e-|?l-o-le-&l-o-l(i7^_^)-ieSlxo-||,-|?|x„-.|^^_ 

The second parenthesis contains a Hilbcrt-Schmidt operator, while the other factors are bounded. 
The decay comes from the last factor. 

Estimates (jB.15|) and (jB.18[) can be proved as the previous ones, taking into account that 
z = —1 can be used in our Combes-Thomas estimates (Proposition [BjT]) , since it belongs to the 
resolvent set of Hd ; see also the proof of Proposition IB. 31 □ 

B.2 Pointwise bounds on the Green's function 

The goal of this subsection is to prove a pointwise bound on the resolvent kernel Q(-,-,z) of 
{Hu — where the operator Hjj is defined as in the previous subsection. 

Proposition B.8. There are constants C > and S > such that for all z Cz V and for all 
X, x' e one has 

|g(x,x',z)| < C-(l + |ln|x-x'||)e-^l"-'''l. (B.20) 
V 

We emphasize that the constants C and S can be chosen independently of A. 

Proof. We use the resolvent identity (|B.10p . By the diamagnetic inequality, as in the proof of 
Corollary El 

|g(x, x', -1)1 < C (1 + I In |x - x'l I) e- /2 . 
Moreover, since (x — x') < (x — y) + (y — x'), 

1/2 



/ dyg(x,y,z)g(y,x',-l) < f / dy e'^<^-y^\g{^,y, z)\A 

/ dye'^(y--'^giy,^\-l)\A 

JR2 / 



1/2 



By (|B.9[) and (|B.7p . both factors are finite and their product is bounded by a constant times 7'^/?^. 
This proves (|B.20p . even with |x — x'| in the exponential replaced by (x — x'). □ 

B.3 Pointwise bounds on the derivative of the Green's function 

In this subsection we specialize to the case of a constant magnetic field and an operator defined 
on all of . We shall prove 

Proposition B.9. There are S,C > such that for any A G C^{M.'^) with curlA constant and 
for any z G T>, the kernel Goi', ■, z) of ((— ?V — A)^ — z)~^ satisfies 

|(-zVx - A(x))eo(x,x', z)| < C-|x - x'r'e-^l^-^'l . 

V 

We emphasize that the constants C and S can be chosen independently of the value of the 
(constant) magnetic field. 
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Proof. By gauge covariance of the kernel we may assume that A(x) = |(— a;2, xi) for some 6 £ M. 
Since the case 6 = is well-known and since the case 6 < can be reduced to the case 6 > via 
complex conjugation, we may also assume that 6 > 0. 

According to Mehler's formula, the heat kernel oi H ^ (— iV — A)-^ is given by 

exp(-ti/)(x,x') - ^exp(-^)exp(-±-^). (B.21) 

Hence 

(-zV. - A(x)) exp(-tH)(x, x') ^ ^ exp(-ii/)(x, x') (^^^^ - (x - x')^ 

where (x — x')^ = (— (^2 — x'2),xi — x'l), and 

|(-iVx-A(x))exp(-<H)(x,x')| = ^^"^ ~ y^tanli~'(bt) + 1 |exp(-tg)(x, x^ 
Using that a > tanh a for all a > and that 

c := sup — - — V tanh^^ a + 1 < oo , 
a>o smha 



we find that 



-^V.- A(x))exp(-tF)(x,x')| < ^^P^^^^^ " 



We will compare this bound with 



X — X 

Vxexp(tA)(x,x') = — — 5-exp(- 



From the last two relations we deduce that for any z with Re z < 



-iVx-A(x))go(x,x',z)| 



/ die*^(-iVx - A(x))exp(-tiJ)(x,x') 
Jo 



oo 

87rt2 At 

-1/ 



<c|x-x'| / dte*'^°^-^exp( 



=c|Vx(-A-Rez)-i(x,x')| . 

Using the explicit form of V(— A — Re z)'"^ as convolution with a Bessel function and the properties 
of that function we deduce that for some C > 



-iVx - A(x))go(x, x', -1)1 < C|x - X' 



/|-lg-|x-x'|/2 



(This is not optimal, but sufficient for our purposes.) 

Now for arbitrary z G 2? we write, using the resolvent identity, 



-iVx - A(x))eo(x, x', z) =(-iVx - A(x))go(x, x', -1) 

+ (z + l) / dyHVx-A(x))go(x,y,-l)go(y,x',z) 

JR2 



This, together with the estimate on ^o(y;x',z) from Proposition [B^ and some simple estimates, 
completes the proof. □ 
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C Existence and continuity of certain integral kernels 

Let H^) be as in the previous section, but defined on the whole space R^. 
Proposition C.l. 

For any f m S(R),, the integral kernel /(i?D)(x,x') of /{Hd) belongs to C"(M.'^ x R^), 
Proof. We start by estabhshing the fohowing lemma. 
Lemma C.2. 

Let U be a compact set in and x be a smooth and compactly supported function which is equal 
to 1 on U. Then, for sufficiently small e > 0, we have that the operator x{'){Hd + l)^^e'^^'^ 
maps _L^(R^) into iJ^(R^). Moreover, for every tp <E L^(R^), the function x{'){Hd + l)^"^e^^'^'0 is 
continuous on U . 

Proof. Let us write: 

x{-){Hd + l)-ie^<-> = [e'^\{-) [Hd + l)"'} [{Hd + l)er'^HHD + l)-^^^^-^} . (C.l) 

A standard argument shows that x{'){Hd + 1)^^ maps L^(R^) into iJ^(R^). The first factor 
e'^^'^ on the left remains bounded because x has compact support. Now if e is small enough, the 
estimate ((R2|) implies that {Hd + l)e~^^'\HD + l)"^e"< > is bounded on L'^{R^). This proves the 
first statement of the lemma. The second statement is a consequence of Sobolev embeddings. □ 

Wc now can start the proof of the proposition. If we introduce f{t) := {t + then we 

write 

x{-)f{HD)x{-) = {x{-){Hd + l)-ie^<->} e-^<->/(iJc)e-^<-> {e^'^HHo + • 

The operator W := e^'^^'^/(iJ£))e~^^'^ is Hilbert-Schmidt, because we can write it as 

W = [e-'< \HD + 1)-^} {Hd + l)/(ifi5)e-'<-> 

and use the fact that c^^^'^Hd + 1)"^ is Hilbert-Schmidt (use ([BJl) ). 

Hence W is an integral operator, with an integral kernel in i^(K"'). 
Therefore we are left with the investigation of the continuity on J7 x [/ of the integral kernel defined 
by: 

/(i/z?)(x,x'):- j j g{y,^;^l)e<^^W{y,y')e<^'^g{y',^-^l)dydy'. 
Using (|B.7p . that W is in and the Cauchy-Schwarz inequality, we obtain 

\f{HD){^,^)\<C\\W\\B,. (C.2) 

uniformly in x, x' € U , where the constant C only depends on the diameter of U . Since W is 
Hilbert-Schmidt, we can approximate W{y, y') (in Hilbert-Schmidt norm) with a finite sum of the 
type J2 ^j(y)^j(y') where the fc's and /I's are i^-functions. Applying the inequality (jC.2[) to the 
difference of the two resulting operators, we see that we can approximate the function f{HD){-, ■) 
(in the uniform norm on U x U) with functions of the type 



Y.{{Hd + l)-ie'<->fc,}(x){(i/z5 + l)-ie^<-)/i,}(x'), 

which from Lemma IC.2I are continuous on compacts. Hence f{HD){-,-) is continuous on its 
variables onU x U. Since U was arbitrary, the proof is over. □ 
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